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ABSTRACT 


This manual describes the use of the computer program BUCLASP3. The 
code is intended for thermal stress and instability analyses of 
structures such as unidirectionally stiffened panels. There are 
two types of instability analyses that can be effected by PAINT; 

(1) thermal buckling and (2) buckling due to a specified inplane 
biaxial loading. Any structure that has a constant cross section 
in one direction, that may be idealized as an assemblage of beam 
elements and laminated flat and curved plate strip-elements can 
be analyzed. The two parallel ends of the panel must be simply 
supported, whereas arbitrary elastic boundary conditions may be 
imposed along any one or both external longitudinal side. 

Any variation in the temperature rise (from ambient) through the 
cross section of a panel is considered in the analyses but it must 
be assumed that in the longitudinal direction the temperature field 
is constant. Load distributions for the externally applied inplane 
biaxial loads are similar in nature to the permissible temperature 
field. 

This manual consists of instructions for the use of the program 
with sample problems, including input and output information. The 
theoretical basis of BUCLASP3 and correlations of calculated results 
with known solutions, are presented in Ref. 1. 
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1.0 SCOPE OF BUCLASP3 

The computer program BUCLASP3 is designed to perform the following types 
of analyses for unidirectional ly-stiffened, rectangular composite 
panels. 

(1) Thermal stress 

(2) Thermal buckling 

(3) Buckling due to a mechanical biaxial inplane load 
distribution varying through the panel cross section. 

The structure must be idealized as an assemblage of beam elements and 
laminated curved and/or flat plate strip elements. Each element 
extends the full length of the structure. The element edges which 
are normal to the longitudinal axis of the panel are assumed to be 
simply supported and any external edge parallel to the longitudinal 
axis may be arbitrarily constrained by specifying appropriate spring 
constants. When needed, plate offsets can be utilized to correctly 
idealize a structure. 

The stress and buckling analyses are based on the linear elastic 
theory. In the buckling formulation, prebuckling deformations and 
initial imperfections are not considered. Applied loads may be 
comprised of axial and transverse plate inplane loads, and N^ 2 j 
and axial beam loads. These inplane applied loads are input as internal 
loads and are assumed to be constant in the longitudinal direction. 

As discussed in Ref. 1, the selection of sinusoidal variations of 
displacements in the longitudinal direction makes it possible, for 
purposes of analysis, to treat the structural elements as if they 
were in the cross section of the structure. Thus, each two- dimensional 
plate element is completely defined as a one-dimensional member and 
each beam element is defined as a point in the. panel structure cross 
section. 
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By prescribing a set of wave numbers, one obtains a separate buckling 
solution for each corresponding half wave length. For each solution, 
the mode shape in the longitudinal direction is sinusoidal and has 
the number of specified half waves. After the solutions for each of 
the input wave numbers are calculated, the minimum buckling load is 
automatically selected and if specified, the mode shape through the 
cross section of the structure is computed for the critical wave 
number. 

In the buckling analyses of (2) and (3) above, the plate elements 
must be subdivided along the width into a specified number of sub- 
elements (see Figure 2,3.1). The division is made such that the 
subelements have equal widths and span the total length of the panel. 

The applied loads are assumed to be constant in each subelement, but 
may vary from one subelement to another. The beam axial loads are 
also assumed to be constant. Shears and moments for all element types 
are ignored in the analysis. 

The thermal stress and thermal instability analyses account for 
temperature variations through the panel cross section. It is assumed, 
however, that the temperature is constant in the longitudinal direction. 
The load distributions for (3) above are similar in nature to the 
permissible temperature field. 



1.1 Thermal Stress Analysis 

Temperature variations within a panel that may be analyzed by the 
thermal stress analysis program were discussed previously. In each 
case, the temperature distribution (a constant value) is expanded as 
a truncated Fourier sine series. By choosing appropriate Fourier 
series expansions as functions of the longitudinal direction for all 
displacements and stresses, the thermal stresses are computed for each 
odd numbered harmonic (wave). The final stresses and displacements 
are aeterminea Dy summing trie contributions from each harmonic. 

These stresses and displacements are printed with respect to element 
local axes. The adopted stress and displacement sign conventions 
for plate elements and beam elements are illustrated in Figure 2.2 
and 2.4.1, respectively. 
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1.2 Thermal Instability Analyses 

Two types of thermal instability analyses can be performed by BUCLASP3. 

(1) Buckling caused by thermal effects only (thermal buckling). 

The critical temperature ratio, defined as the factor by 
which the specified temperature distribution is multiplied 
where buckling occurs, is computed. 

(2) Buckling due to mechanical axial loads, (mech.), on a 
panel exposed to a specified invariant temperature distribution. 

In both cases, a thermal stress analysis is performed first. In general, 
the computed thermal plate stresses, (therm.) and (therm.) are 
characterized by a two-dimensional variation within each subdivided 
plate element. These stresses are calculated at the subelement mid- 
widths and are averaged in the longitudinal direction. These average 
stresses are then applied to each subelement as constant prestress loads. 
All other thermal stress components are ignored in the buckling analysis. 
The prestress loads, constant within any subelement, but varying from 
one subelement to another, account for the possible varying nature of 
the plate thermal stresses in the transverse direction. Beam element 
axial loads are averaged and are applied as constant prestress loads 
for the buckling analysis. 

In (1) above, the factor that all plate subelement and beam prestress 
loads are multiplied by, which causes buckling, is computed. This 
factor, greater than 1.0, is called the critical temperature ratio. 

At buckling, the panel prestress load distribution is identical, with 
the exception of magnitude, to the initially computer prestress load 
distribution. 
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For (2) above, a mechanical axial load, (mech.) is superimposed 
on the prestress loads (from thermal stresses). The (mech 0 ) axial 
load is distributed to all elements such that longitudinal strain 
compatibility is satisfied. This compatibility condition is satisfied 
with the (mech.) load only, excluding the prestress loads. The 

critical (mech.) load that causes panel buckling is computed for 

this option. The thermal prestress loads remain invariant throughout 
the buckling computations. 

For both options above, the user must ensure that the input temperature 
distribution is not too large that it hasn't already caused the panel 
to buckle. This is necessitated by the restrictions imposed on; 

(1) above, where only critical temperature ratios greater than one 
are considered, and (2) above, where tensile (mech.) loads are 
excluded from the calculations. 


1.5 



1,3 Instability Analyses of Panels Subjected to Mechanical Loads 

There are two options for buckling analysis of panels subjected to 
mechanical, biaxial inplane loads. They are: 

(1) Critical load ratio option, where the critical load ratio 

is defined as the magnitude of the buckling load distribution 
divided by the magnitude of the initially specified inplane 
load distribution. 

(2) Buckling option due to an initially specified load distribution 
upon which an additional applied axial load, is super- 
imposed. The critical value of is computed for a panel 
that is initially prestressed by input loads. The prestress 
loads are assumed to remain constant throughout the buckling 
computations. 

These two options are treated in an identical manner as the options 
previously discussed for thermal instability analyses. Replacing the 
computed thermal prestress loads by the input plate subelement loads, 

N-|i (input) and (input), and the input beam axial loads, the 
buckling analyses are performed in the same manner as the thermal 
instability analyses. 

Again, only uniform N.^ and N ^2 inplane loads for plate subelements and 
constant axial loads for beams are considered. The same restrictions 
defined for the thermal instability options are also applicable here. 
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2.0 BASIC INFORMATION 

This section presents basic information required by BUCLASP3 users 
for performing stress and buckling analyses. Figure 2.1 shows 
the cross section of a structure. The circled numbers are element 
numbers. Extremities of each plate element and the shear centers 
of each beam element are indicated by dots. The dashed line is 
drawn through the mid-planes of each plate element. 
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2.1 Geometry 


2.1.1 Coordinate Systems 

All structural geometries are defined relative a global (overall) 
right-handed rectangular Cartesian coordinate system. The axes 
of this triad are denoted by Xg, Yq and Zq in Figure 2.1. 

The longitudinal axis of the panel is in the same direction 
as X fi . Local element coordinate systems will be discussed 
individually for each element in Section 2.2. 

2.1.2 Nodes 

The dots shown in the Figure 2.1 will be referred to as nodes. 
They are actually lines that are parallel to Xq . A node is 
identified by its coordinates in the global system and by its 
user-assigned node number. This node number is only for the 
user's convenience since the program automatically numbers the 
nodes sequentially according to the input order. The first input 
node will be internal node number 1, the second input node, 2, 
etc. Generally, the user-assigned node numbers will not be 
the same as the internal node number. No reordering of nodes 
is done internally. 

2.2 Elements 

There are three basic types of elements: (1) flat plate, (2) 

curved plate, and (3) beam elements. All elements extend the 
full length of the structure. A plate element is defined by a 
pair of node numbers that correspond to the sides of the element. 

In Figure 2.2 the plate element is defined by nodes A and B. A 
beam element is defined by a node that corresponds to the beam 
shear center. In Figure 2.1 the following observations can be made 
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(a) CD , ® , © , © . ® . and © are flat plate elements. 

(b) (2) , ® , and ® are curved plate elements. 

(c) © and © are beam elements. 
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2.2.1 Plate Elements 

There are two types of plate elements, flat and curved. The 
following discussion applies to both types except when reference 
is made to a particular type of plate element: 

Local Coordinate System 

The local coordinate system of a plate element is defined by 
the ordering of the node numbers that identify the element. 

The local y axis of the plate element is in the direction from 
the first node of the element to its second node. For example, 
in Figure 2.2 the element is defined from node A to node B, and 
the resulting direction of the local y axis of the element is 
as shown. The local y axis of the curved plate element is along 
the arc from A to B. The local x axis is parallel to and in the 
same direction as Xq. The local z axis is determined by the 
right-hand rule. The origin of the local axes is on the mid- 
surface of the element and midway between its nodes. 

Properti es 

In its most general form, the plate element can be multi-layered, 
each layer being orthotropic with respect to the local plate 
axes. In its simplest form, the element degenerates to a single 
isotropic layer. The thickness of each layer is assumed to be 
constant. The user must specify the thickness and material 
property of each layer. 

The first lamina of a laminated plate element is defined as the 
external layer in the negative z direction. The last lamina is 
the extreme layer in the positive z direction (see Figure 2.2). 

It is important that the layers are input in this particular order. 
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Radius of Curvature of Curved Plate 

The curved plate element is actually a strip of a circular 
cylindrical shell. Therefore, only one radius of curvature, 

R, need be defined for this element. The sign convention for R 

is defined in the following manner: R is positive if the vector 

from the element to the center of curvature is in the same 
direction as the local z axis of the element. The magnitude of 
R is measured from the mid-surface of the element to the center 
of curvature. In Figure 2.3 the sign of R for elements (a) and 
(b) is negative, whereas it is positive for element (c). 

Temperature Distribution 

The temperature field within a plate element is expressed as: 

T p = T px' T py' T pz 

where Tp = plate temperature, 

T px = variation of the temperature field in the x direction, 

T = variation of the temperature field in the local y 

py direction, 

and T = variation of the temperature field in the local z 
p direction. 

Tp x is assumed to be a constant as discussed previously. Tp^ may 
be specified either as a constant or as a variable by inputting 
the coefficients, a^, b^ and c^, in’ the following polynomial, 

T py = V 2 + V * V 
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where y is the local axis of the element in the plate width 

direction. 1. , which is treated in the analysis as a constant 
pz 

for each layer of the plate element, may be specified in two 
different manners; (1) a constant value for each layer or (2) 
a variable expressed by the coefficients A z> B z and C z in the 
following polynomial. 



V + C z z 
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In this relation, z is the distance measured along the local z 
axis of the plate element from the element outer surface. When 
a user defines Tp Z as a function of z, the program computes Tp Z 
at the mid-surface of each layer and applies this value to the 
layer as a constant. 


If the thermal field varies in the local z direction within a 
layer, the layer may be idealized in the following manner. 
Assume this lamina to be multi-layered and input constant 
values for Tp Z to approximate the actual variation through the 
thickness. 


Applied Loads 

When the option to compute buckling loads for a panel subjected 
to only mechanical loads is used, the initial inplane biaxial 
loads must be specified. They are specified as constant values 
for each subelement. Tensile loads are input as positive values 
and compression loads as negative values. 


2.9 




Figure 2.3.1 Subdivided Plate Elements 
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2,2.2 Beam Elements 

There are two basic types of beam elements that are available 
in BUCLASP3. They are the laminated beam element, rectangular or 
circular in cross section, and the general beam element. The 
local x axes of all beam elements are directed in the same 
direction as the local Xq axis and parallel to it. 

Laminated Beams 
Rectangular: 

The rectangular laminated beam is shown in Figure 2.4(a). The 
depth of each layer is b. The origin of the local axis system 
is at the shear center, 0, of the beam. The local y axis of this 
element could either be in the direction as shown in Figure 2.4(a) 
or in the opposite direction. The local z axis is determined by 
using the right hand rule. 

The laminae are ordered sequentially in the direction of the local 
y axis as indicated by Figure 2.4(a). Thus, if the local y 
axis were in the opposite direction, the ordering of the layers 
would be reversed. 

The temperature distribution within each beam layer is expressed 
as: 


T b = T bx‘ T by' T bz 
where T^ = beam temperature, 

T. = temperature variation in the axial direction, 

T^ = temperature variation in the local y direction 
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C - Neutral Axis 
0 - Shear Center 


(a) Laminated Rectangular Beam 


(b) Laminated Circular Beam 


Figure 2.4 Laminated Beam Elements 
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and T bz = temperature variation in the local z direction. 


T bx and are assumed to be constant and T bz is specified 
either as a constant for each layer or by the following polynomial. 


T bz = a z + b z z + V 


where z is the local axis of the beam. The T bz values are 
automatically calculated at the middle of each layer and these 
values are used in the analysis as constant values. 


Circular: 

The beam is laminated in the manner shown in Figure 2.4(b). 

The first layer is the innermost layer and the ordering of the 
layers extends outward. The local y axis of this element can 
be in any direction in the plane of the beam cross section. 

The local z axis is again determined by the right hand rule. 

The temperature distribution is specified as a constant for each 
lamina. 


General Beam Element: 

The local y axis of the element can be in any direction in the 
place of the beam cross section. The local z axis is defined 
by the right hand rule. Properties such as inertias, areas, 
torsion constants and warping constants are input with respect to 
the local coordinate system. The orientation of a beam element 
is defined relative to the global coordinate system by specifying 
the angle between the global and local axes. 

For the thermal analyses, the thermal loads are input for this 
beam type. These thermal loads can be computed from equations 
(4.53) to (4.55) of Reference 1. 
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2.3 Offsets 


The element stiffnesses are originally defined at the element "local 
nodes". A plate element has two "local nodes" defined as the longitu- 
dinal edges lying in the mid-surface and the "local node" for a beam 
element is defined by its shear center. When two or more elements 
are joined together, the appropriate element stiffnesses can be 
merged at the juncture if the "local nodes" of all elements are 
coincident at the juncture. If the "local nodes" at this juncture 
are not coincident, it is necessary to define a common node where all 
the element stiffnesses can be merged properly. By necessity this 
common node must be a user input node. By using plate offsets, the 
user effectively defines the element stiffness at some node other 
than the plate "local node", thereby, giving him a tool to meet the 
above requirements. This non-" local node" is referred to as an 
"offset node". It must, however, be one of the user-defined nodes. 

A plate with offsets can be interpreted as a pseudo plate element with a 
"link" connecting one of the plate "local nodes" to an "offset node". 

This link allows the element stiffness at a "local node" to be transferred 
to an "offset node". 


There is no offset capability for the beam element in this program. This, 
however, is not detrimental since in all cases plate offsets can be used 
to idealize non-coincident junctures. 

A situation where plate offsets can be used is when the mid-surfaces of 
plates do not intersect at the juncture. Offsets can also be used when 
the shear center of a beam does not coincide with the mid-surface of 
an adjoining plate. 
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The plate offset is a vector which is specified by two components, y Q 
and z Q from the edge of the plate element (see Figures 2.5.1 and 2.5.2). 
The y Q and z q vector components are measured relative to the local 
y and z axis respectively. They are positive when their directions 
coincide with the positive directions of the local axes, z is 
measured from the negative-z plate surface (exterior surface in the 
negative-z direction) and y is measured from the plate "local node". 

It must be emphasized that z q is not measured from the "local node". 

Figure 2.5.1 illustrates the cross section of a panel in which two plate 
elements intersect a beam element. The input nodes coincide with the 
plate "local nodes" but they are not coincident with the beam shear 
center. Since these nodes are not coincident at the juncture and a beam 
offset capability is not available, offsets for plate elements (T) and 
(T) must be specified. Consequently, the element stiffnesses will be 
merged at 3, the beam node. 

Figure 2.5.2 illustrates two examples where plate elements of unequal 
thicknesses meet. In the top figure, the mid-surfaces of plate elements 
(T) and (3) with equal thicknesses passes through the input nodes. 

Element (T) , however, is a two-layered plate with a different thickness. 
Since the stiffness for element ( 2 ) is formulated relative to its "local 
nodes", offsets must be defined to transfer this stiffness to nodes 2 and 
3 to ensure the desired structural integrity. If offsets were not 
specified for this example, the true configuration of this structure 
would not be retained in the analysis. That is, element (J)would be 
situated such that its mid-surface coincides with the mid-surfaces of 

©and(3> 

It should be noted that z is not the offset distance from a plate 
"local node" to the "offset node". This offset component, however, is 
used by the program to formulate the actual offset-transfer of the 
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stiffness. For example, the stiffness at the mid-surface of (J) in 
Figure 2.5.2, is transferred to nodes 2 and 3 based on the z q offset 
indicated therein. 
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y Q Positive ( for Element (T) ) 
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Figure 2.5.1 Offsets in Presence of Beams 
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Z.4 Boundary Conditions 

An external side parallel to the Xq axis of a plate element may 
be constrained by specifying one of the following types boundary 
conditions : 

(a) .Simply supported 

w = M «2 = ^22 = u = 0 

(b) Clamped 

w = 0 = N 22 = u = 0 

(c) . Free 

n = M 2 2 = N 22 = N 12 = 0 

(d) Sprung 

The constrained or sprung degrees of freedoms are defined with 
respect to the local coordinate system of the plate element. 

Boundary conditions cannot be applied to a beam element node or 
an internal node. For a sprung node of a plate element, the 
spring stiffnesses for all four degrees of freedom must be specified. 
To delete or constrain a freedom, use a large value for the 
corresponding spring constant, e.g., 1 X 10^. 
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2.5 Substructures 

Three types of substructures are considered in PAINT; (I) start, 

(.2) repeat and (3) end. The basic reason for including these 
substructures is to avoid repetitious computations that is 
generated by panels with many repeated structural groups. Also, 
higher efficiency can be achieved by dividing the stiffness matrix 
into separate blocks, each block corresponding to a substructure, 
during the buckling load computation. The user input is also 
reduced. 

The user may divide a panel, with or without repetitive structural 
groups, such that is consists of one start substructure, one, many, 
or no repeat substructures and one end substructure. Further, a 
panel may be idealized with the following combinations of substruc- 
tures; 

(1) start substructure only 

(2) start and end substructures only 

(3) start, repeat (one or more) and end substructures. 

It should be obvious at this point that panels with no repeatable 
structural groups can have a maximum of one repeat substructure. 

It is recommended that the start substructure be larger of the 
three substructures, for efficiency reasons in the determinant 
cal cul ati ons . 
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2.5.1 Flat Panel with Repeated Substructures 

For a panel with, repeated substructures, a reduced version of the 
actual panel may be specified for input. For example, see Figure 
2.7.. (a) is the actual panel and (b) is the reduced panel used 
for input. The nodal coordinates of 7 to 10 of (b) are input such 
that they model the end substructure of the actual panel, (a), 
as if it is moved to the left and attached to the first repeat 
substructure. Therefore, nodes 7 to 10 of the reduced panel, (b), 
are related to nodes 23 to 26 of (a), the actual panel. Consequently, 
nodes 11 to 22 of (a) are not directly included by coordinates in 
the input, but are implicitly specified through the repeat sub- 
structure. 

After drawing a sketch of the actual structure, the user must 
identify the nodes that belong in the three different substructure 
types. For the start and end substructures, the only limitation is 
that there must be at least one interior node, a node without any 
specified boundary conditions, in each substructure. To ensure 
that the repeat substructure is accurately identified, the sub- 
structure corresponding nodes must be identified. The corresponding 
nodes are those node pairs at which the stiffnesses, independent 
of element temperatures and loads, are identical . In other words, 
all elements meeting at corresponding nodes in all repeat sub- 
structures must be identical in all respects, including boundary 
conditions, except for the element temperatures and element loads. 

For the example from Figure 2.7(a), inspection shows that nodes 7, 

11, 15 and 19 are the corresponding nodes to node 3. Therefore, 
note 3 is a valid node for the repeat substructure. Similarly, 
node 6 is a valid node for the repeat substructure with nodes 10, 

14, 18, 22 as correspondence nodes; node 4 with 8, 12, 16, 20; and 
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node 5 with 9, 13, 17 and 21. One should note that node 23 is 
not a correspond! nq node of node 3 because the specified boundary 
condition at node 26 results in the stiffness contribution of 
element 23-26 to node 23 being unequal to the stiffness contribution 
of element 3-6 to node 3. Thus node 23 may not be included in 
the repeat substructure. 

From this user exercise one can conclude that, a) the substructure 
unit made up of nodes 3, 4, 5, 6 is a repeat substructure, and b) 
there are (5) repeat substructures. 

Node 2 is included in the start substructure because boundary 
conditions are specified for node 1. The selection of nodes for 
the end substructure is obvious at this point. 

Although only a reduced version of the actual panel is input, the 
temperatures or input prestress loads for every element of the 
actual panel must be specified. This includes all elements in all 
reoeat substructures. At this point, all elements must be assigned 
to different substructures. Any plate element with both definition 
nodes lying within a substructure is assigned to that substructure. 
If one node of an element is in the start substructure and the other 
node in the repeat substructure, the element is included in the 
repeat substructure. Plate elements, with one node in the repeat 
substructure and the other node in the end substructure, are 
assigned to the end substructure. Since a beam element is defined 
by a node, the assignments of beam elements should be obvious. 

In figure 2.7(b), element (7) is in the start substructure, elements 
© » (J) » ©) cmd(© are in the repeat substructure and elements 
© . © . ® and © are in the end substructure. 
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2 0 5 e 2 Substructure Interrelationship Plate Element Pairs 

The mode shapes for elements of the panel (not the reduced 
structure) is provided by the program. These elemental mode 
shapes are derived from the relative displacements which were 
determined individually for each node of the actual panel 
(not just the reduced structure). The previous input information 
on repeat substructures was adequate to produce (internally) 
the total structure stiffness from the (input) reduced structure. 
However, additional information is required to identify elements 
that connect one repeat substructure to another repeat substructure 
for the purpose of calculating elemental mode shapes. 

From the definition of a repeat substructure , any element, that 
connects the start substructure to the adjoining repeat substructure, 
has a corresponding element (topologically similar), that connects 
the repeat substructure to the end substructure. The requirement 
then becomes one of identifying; 

(a) Each connecting element between the start and repeat 
substructures 

and 

(b) each connecting element (topologically corresponding to (a)) 
between the repeat and end substructures* 

These elements become logically an element pair, and are defined 
as "interrelationship element pairs". 

As an aid to the user, the following steps are recommended (see 
Figure 2.8). 

(1) After encircling the nodes for different substructures 
as shown in Figure 2*8, identify all element(s) that 
connects the start substructure to the repeat substructure. 
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(Elements (T) , (?) , (§) , (6) and (7) of Figure 2.8) 

(2) Identify all element(s) connecting the repeat and end 
substructures. (Elements © , © , © , © and 
© of Figure 2.8) 

(3) Pair all elements of (1) above, with corresponding elements 
(topologically similar) from (2) above. For the panel in 

Figure 2.8, the element pairs are: ( (3) , © ), ( © > © ), 
( © , © )» ( © > © ) and ( © , © ). These are the 
interrelationship element pairs. 

As an additional example, see Figure 2.7(b). The interrelationship 
element pair is ( © , © ). 
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Figure 2.8 Example for Repeat Substructure Interrelationship Element Pairs 
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2„5.3 Some Restrictions on Substructuring Non-Repeti ti ve Panels 

Whenever a panel without repetitive structural groups is 
idealized with more than one substructure type, care must 
be exercised in the substructuring process. The panel must 
be divided such that any one substructure contains only one 
unconnected piece of structure. To illustrate this, Figure 2.9 
is presented. This skin-stringer assembly is divided into 
three substructures. Node 3 is not attached to node 4, node 6 
is not connected to node 7 and node 5 is used to attach the 
stringer to the skin. In Figure 2.9 (a), the repeat and end 
substructures defy the restriction stated above, because in 
the repeat and end substructures, there are two discontinuous 
or unconnected pieces of structures. In the end substructure, 
the structure containing nodes 9 and 11 is not attached to the 
piece of structure that contains nodes 12 and 10, and in the 
repeat substructure, the piece of structure containing nodes 7 
and 8 is not attached to node 6. Therefore, there are two 
unconnected structural pieces in each of the two substructures, 
repeat and end. Figure 2.9 (b) is an example of a valid 
selection of nodes for this problem. Each substructure contains 
only one continuous piece of structure. 

Another restriction is in the input order of the nodes. The 
nodes should be ordered such that the input list of nodes can 
be separated into distinct groups, each group related to a 
different substructure. The start substructure nodes should 
be input first, repeat substructure nodes, second, ard the end 
substructure nodes last. What is important here is the input 
or order of nodes, not the user node numbers. If the nodes of 
Figure 2.9 (b) are input in the same order as the user node 
numbers, an incorrect data set will result. Nodes 9 through 11 
would be out of order. Instead, the input order should be 
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(1,2,3,4,5,6,7,8,9,11,12,10). The first group of nodes for 
the start substructure would be (1 to 6) , second group for the 
repeat substructure , (7,8,9,11) and the last group for the end 
substructure, (12 and 10). 
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Figure 2. 9 Examples of Node Restrictions in Substructures 
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2.6 Upper and Lower Bound Loads 

The panel upper bound load is required to increase the efficiency 
of the method used in extracting the lowest eigenvalue. This 
panel upper bound is defined as the least of the element upper 
bound loads. An element upper bound load is the buckling load 
of a plate element with interior sides completely restrained. 

For a plate element with specified boundary conditions on one 
side, the upper bound load is the buckling load for the plate 
with the interior side completely restrained and the exterior 
side subjected to the specified boundary conditions,. 

The upper and lower bound loads can be specified as input to 
PAINT. If the bounds are not specified, an approximate upper 
bound load will be computed. Since all plate elements are subdivided, 
the subelement upper bound loads must first be determined. For an 
interior subelement, they correspond to buckling loads for a subelement 
completely restrained at both sides. These loads are computed by 
using the Galerkin method as described in Ref. 2. For an exterior 
subelement (subelement with a side with specified boundary conditions), 
the buckling determinant is first formed and then the buckling load 
is solved. The minimum of the subelement loads is selected to be the 
upper bound for the element. Finally, the element upper bound is 
calculated by setting the element buckling determinant to zero. The 
initial prestress loads are chosen as the element lower bound. 
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2.7 Mode Shapes 

The relative displacements may be computed for the buckling 
load corresponding to the critical wave number. The relative 
displacements are computed for each element with respect to 
its local coordinate system. They can be computed at evenly 
spaced prescribed points along the width of any plate element. 
From the output of relative displacements, the mode shape can 
be hand plotted. 
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2.8 Modelling of Riveted Connections 

This section deals with the modelling of riveted connections. 

As an example, see Figure 2.10. If the rivet is ignored in 
the idealization, the stiffener would be assumed to be bonded 
to the skin. A suggested modelling technique is presented 
here to more correctly depict the behavior of riveted connections. 

The rivet is idealized as a line at nodes 4 and 5. In Figure 2„10(a) 
the plate elements (J) and (T) are offset at node 4 as shown. Node 
5 is equivalenced to node 4 to create an attachment between these 
nodes. The plate elements (2) , ( 3 ) , ( 7 ) , and(5) are all joined 
at node 4, the rivet line, but node 2 of element (?) and node 3 of 
element (T) are not connected and they are free to deform independently. 
Nodes 6 and 7 of elements ( 4 ) and (s) respectively, are also free to 
move independent of each other. This modelling is approximate in 
the sense that the connection between the stiffener and skin is 
reduced to a line connection. 

First, make an initial run with this idealization shown in Figure 
2.10(a) and investigate the mode shapes to see if they are feasible. 

For an example of mode shapes not being feasible, assume that the 
downward deflection of node 2 in the direction is greater than 
the downward deflection of node 3 in the same direction. Wherever 
this type of physical unreality occurs, join the elements together at 
these nodes and rerun the problem. For the example cited above, the 
connection of nodes 2 and 3 is accomplished by making node 2 equivalent 
to node 3. Offset element(3)as indicated in Figure 2.10(b). 
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In the original idealization shown in Figure 2.10(a), nodes 
2 and 6 seem unnecessary for the first trial analysis. Plate 
elements (T) and (T) and similarly elements (4) and (8) could 
be combined into one element. The only reason for including 
nodes 2 and 6 in the example is to compare the displacements 
at these locations with nodes 3 and 7 of the stiffener, to see 
if any unrealistic deflections occur at these points. Actually, 
the displacements along the junction line between the skin and 
stiffener top flange should be investigated to check the feasibility 
of the deformed shapes. 
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3.0 PROGRAM COMPUTER DETAILS 


3.1 Machine Requirements 

The BUCLASP3 program is written for the CDC 6600 series computers. It 
requires the use of a card reader, line printer, disk storage, a 
minimum of zero and a maximum of one tape drive. 

3.2 Operating System 

The program runs under the SCOPE 3.1 or KRONOS 2.0 operating systems. 
All system routines used are assumed to be standard CDC release. 

With the exception of four special purpose routines in COMPASS, all 
source routines are coded in CDC FORTRAN IV. The overlay loading 
feature is used. 

Two Langley Research Center routines are used. 

JPARAMS COMPUTER PROGRAMMING MANUAL 

Volume I, Section K5.1 date 6-9-71 

This routine makes certain job parameters available to the user. 

SIMEQ COMPUTER PROGRAMMING MANUAL 

Volume I, Section F4.1 date 8-1-68 

This routine solves a set of linear equations. 
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3.3 Storage Allocation 

The program will LOAD and EXECUTE with a field length of 77000g. 

The field length estimate for the program is dependent upon the 
data set. The space for the buckling determinant is dependent upon 
the size of the three blocks of the determinant definition. 

When the eigenvector and relative displacements or thermal stresses 
are also wanted the user must take into account in field length allo- 
cation that the start substructure, each repeat substructure and the end 
substructure have to be stored in core. 

The data processing overlay prints an estimate of field length required 
for the other overlays. 


3.4 Timing and Output Estimates 

Time consumption for one data set depends on various factors. 

a. Number of wave numbers investigated. 

b. Number of nodes, plate and beam element defined for panel. The 
size of buckling determinant is proportional to number of uncon- 
strained nodes. 

c. Number of plate elements with stiffnesses which are the same. 

d. Number of Fourier Harmonics in thermal stress analysis. 

e. Number of subelements. 

In general, for single data set runs, and with the intermediate 
printout switches off, the program will generally generate from 30 to 
100 pages. The number of wave numbers is a primary factor on the 
amount of printout. The intermediate printout should be used with 
care as it can be voluminous. 
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3.5 Disk File and Tape Utilization 

BUCLASP3 uses a maximum of nine files in addition to the standard input 
and output. The nine files are scratch files used for internal 
information transfer. 


The files referenced are: 


(1) 

INPUT 

Standard Input (Card Reader) 

(2) 

OUTPUT 

Standard Output (Line Printer) 

(3) 

TAPE1 

Scratch 

File 

(Binary) 

(4) 

TAPE2 

Scratch 

File 

(Binary) 

(5) 

TAPE 1 1 

Scratch 

File 

(Binary) 

(6) 

TAPE12 

Scracth 

File 

(Binary) 

(7) 

TAPE! 3 

Scratch 

File 

(Bi nary) 

(8) 

TAPE 1 4 

Scratch 

File 

(Bi nary) 

(9) 

TAPE! 5 

Scratch 

File 

(Binary) 

(10) 

TAPE16 

Scratch 

File 

(Binary) 

OD 

TAPE17 

Scratch 

File (Binary) 


3.6 Control Cards 

There are basically three modes of execution, from 

(a) source 

(b) relocatable binary 

(c) absolute binary 

In the following, use of specific control cards has been avoided, 
rather the required sequence to operations is specified. All file 
names with the exception of BUCLASP are arbitrary. All overlays have 
the name BUCLASP, thus a file BUCLASP is generated at load time. For 
the cases above 
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(a) 1. Obtain a source file, PROG, from permanent storage 

(Cards, TAPE, permanent disk file, etc.). 

2. Compile source file placing relocatable binary on BPRG. 

3. Load BPRG. 

4. Execute from BUCLASP. 

(b) 1. Obtain a relocatable binary file, BPRG, from permanent 

storage. 

2. Load BPRG. 

3. Execute from BUCLASP. 

(c) 1. Obtain an absolute binary file, BUCLASP, from permanent 

storage. 

2. Execute from BUCLASP. 
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4.0 PROGRAMMED DIAGNOSTIC MESSAGES 

There are two types of error messages. The first is a message followed 
by a routine name and number. This kind is fatal to further processing 
so an exit is taken. The second kind is used in the input data processing 
routines to indicate a data error has been detected, but an attempt will 
be made to process the rest of the input data. This kind has a message 
but no number. 


Routines which print diagnostic messages: 


DATAPRO 

FINDIN 

BLKDEF 

TSTIFF 

SMOC 

SOLVE 

TDISPL 

TSTRESS 

LCNTRL 

UPRBND 

DM 

GAPUPP 

DB 

SBLKDT 

DBLERT 

STORE 

PRDMTX 

PROOT 

SOLVV 

DISPLAC 

FETCHD 

READTR 

FNLDIS 


Input Data Processing Routines 
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DATA PRO 


802 

TOO MANY 

NODES 


Card 4 

Field 1 

804 

TOO MANY 

ELEMENTS 


Card 4 

Fields 2,3 and 

806 

TOO MANY 

BEAMS 


Card 4 

Field 2 

808 

TOO MANY 

PLATES 


Card 4 

Fields 3 and 4 


810 BAD EQUIVALENT NODE DATA 

E Card 

812 BAD INTERRELATIONSHIP ELEMENT DATA 

C Card 

901 NUMBER OF VALUES, III, IN M-LIST EXCEEDS 30 
Cards 6 and 7 

902 END OF FILE ENCOUNTERED BEFORE L 
L Card was not found 

903 NON ELEMENT DATA WAS ENCOUNTERED OUT OF ORDER 

904 NUMBER OF NODAL SPRING STIFFNESS SETS EXCEEDED NUMBER 
OF NODES IN DATA SET 

905 NUMBER OF PLATE ELEMENT OFFSET CARDS EXCEEDED THE 

NUMBER OF PLATE ELEMENTS IN DATA SET , 


4.2 



906 NUMBER OF REPEAT SUBSTRUCTURE INTERRELATIONSHIP ELEMENT 
PAIRS INPUT EXCEEDED MAXIMUM OK IIIII 

907 NUMBER OF EQUIVALENT NODE DEFINITION PAIRS INPUT 
EXCEEDED MAXIMUM OF IIIII 

903 INPUT MATERIAL PROPERTY TABLE FAILED 

909 INPUT OF THICKNESS TABLE FAILED 

910 NUMBER OF NODE DEFINITIONS INPUT, XXXXX, DID NOT AGREE 
WITH NUMBER SPECIFIED IIIII 

911 NUM3ER OF ERRORS DETECTED IN PANEL SUBSTRUCTURE 
DEFINITION = IIIII 

912 NUMBER OF PLATE ELEMENTS INPUT EXCEEDED NUMBER SPECIFIED 

914 INVALID LOAD OPTION, IIIII, SPECIFIED 
Card 5 Field 1 

915 INVALID WAVE NUMBER SEARCH OPTION, IIIII, SPECIFIED 
Card 5 Field 2 

916 IIII ERROR(S) DETECTED IN DATA 

917 CURVED PANEL REPEAT SUBSTRUCTURE INTERIOR ANGLE EXCEEDS 
180 DEGREES 

Card 4 Field 6 

918 CURVED PANEL REPEAT SUBSTRUCTURE INTERIOR ANGLE MULTIPLIED 
BY NUMBER OF REPEAT SUBSTRUCTURES EXCEEDS 360 DEGREES 

Card 4 Field 6 
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019 MATERIAL TABLE INPUT, BUT NO THICKNESS TABLE PRESENT 

THICKNESS TABLE INPUT, BUT NO MATERIAL TABLES PRESENT 

NUMBER OF MATERIAL TABLES INPUT EXCEEDS MAXIMUM OF III 

INVALID PLATE TYPE 
Card P Field 5 

NUMBER OF LAM I NAS EXCEEDS MAXIMUM OF 25 
Number of Plates Laminas 25. 

Card P Field 7 

INVALID PLATE OFFSET SWITCH, III 
Card P Field 6 

INVALID MATERIAL PROPERTIES INPUT FORMAT OPTION, III 
Card P Field 8, or 

Card B Field 6 

INVALID N 22 LOAD SELECTION SWITCH, III 
Card P Field 9 

NUMBER OF RELATIVE DISPLACEMENT SUBDIVISIONS WAS NEGATIVE 
Card P Field 10 

INVALID MATERIAL TABLE INPUT OPTION, III 
Card P Field 11 
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INVALID BEAM TYPE 
Card B Field 4 

NUMBER OF LAYERS EXCEEDED MAXIMUM OF 35 
Number of Beam Layers £ 35 

NUMBER OF NODAL CARDS READ EXCEEDS NUMBER DECLARED 

ILLEGAL BOUNDARY CONDITION, 1 1 II I ON NODE III 

NUMBER OF S CARDS, I III I, DOES NOT AGREE WITH NUMBERS 
OF SPRUNG NODES II I II 

USER ID, IIIII, ON S CARD, III II , INVALID 

USER ID, IIIII, IN REPEAT SUBSTRUCTURE INTERRELATIONSHIP 
PAIR, IIIII, IS INVALID 

USER ID, IIIII, IN EQUIVALENT NODE PAIR, IIIII, IS INVALID 
Card E 

USER ID, IIIII, IN, AAAAA, SUBSTRUCTURE DEFINITION IS INVALID 
Card D 

NODE1 ( I I I I ) DOES NOT PRECEDE NODE2(IIII) OF, AAAAA, 
SUBSTRUCTURE DEFINITION 

Card D 

END NODE ( 1 1 1 1 ) OF, AAAAA, SUBSTRUCTURE IS NOT CORRECTLY 
RELATED TO START NODE(IIII) OF, AAAAA, SUBSTRUCTURE 
Card D 
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NODE ( 1 1 1 1 ) FOR PLATE(IIII) IS INVALID 
Card P Fields 3 or 4 


NODE(im) FOR BEAM ( II 1 1 ) IS INVALID 
Card B Field 3 

USER ID, II1II, FOR PLATE ELEMENT OFFSET SET, III II, 

IS INVALID 

Card F Field 2 

NUMBER OF PLATE ELEMENTS INPUT (IIIII) DOES NOT AGREE 
WITH NUMBER SPECIFIED 

NUMBER OF BEAM ELEMENTS INPUT (IIIII) DOES NOT AGREE 
WITH NUMBER SPECIFIED 

minim, is invalid node number 

NODE I, I II I III I II, IS NOT IN START SUBSTRUCTURE 

NODE J, III I III III, IS NOT IN REPEAT SUBSTRUCTURE 

NODE, IIIII, IN START SUBSTRUCTURE SHOULD BE PART 
OF A REPEAT SUBSTRUCTURE INTERRELATIONSHIP PAIR 
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BLKPEF 

708 mil, ERRORS DETECTED IN SUBSTRUCTURE DEFINITION NODE 
PAIRS 

710 PANEL WITH REPEAT SUBSTRUCTURE BUT NO END SUBSTRUCTURE 
DEFINED 

711 SUBSTRUCTURE DEFINITION PAIR IS OUT OF ORDER 

713 NO SUBSTRUCTURE DEFINITION NODE INPUT 

714 PANEL MATRIX DEFINED IS NOT SQUARE 

715 START COLUMN DEFINITION OF SUBSTRUCTURE AA.AA.AA AND AAAAAA 
SUBSTRUCTURE IS NOT RELATED PROPERLY 

Examine Substructure Definition and Element Connectivity. 
See Section 2.6.4. 

716 END COLUMN DEFINITION OF SUBSTRUCTURE, AAAAAA, AND AAAAAA 
SUBSTRUCTURE IS NOT RELATED PROPERLY 

Examine Substructure Definition and Element Connectivity. 
See Section 2.6.4. 

717 ROW DEFINITION OF SUBSTRUCTURE, AAAAAA AND AAAAAA, 
SUBSTRUCTURE IS NOT RELATED PROPERLY 

Examine Substructure Definition and Element Connectivity. 
See Section 2.6.4. 

718 SUBSTRUCTURE AAAAAA AND SUBSTRUCTURE AAAAAA ARE NOT 
INTERRELATED PROPERLY 

Examine Substructure Defintion and Element Connectivity. 
See Section 2.6.4. 
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TSTIFF 

1 


SMOD 

5 

SOLVE 

8 

TDISPL 

1 


2 

TSTRESS 

1 


INSUFFICIENT BLANK COMMON AREA IS AVAILABLE TO CONSTRUCT THE 

MERGED STIFFNESS MATRIX. INCREASE FIELD LENGTH BY 

(OCTAL). 


SINGULARITY DETECTED FOR STIFFNESS SUBMATRIX BEING REDUCED 


SINGULARITY ENCOUNTERED IN ATTEMPT TO PRODUCE STIFFNESS MATRIX 


INSUFFICIENT BLANK COMMON AREA IS AVAILABLE TO CONSTRUCT THE 
COMPACTED MERGED STIFFNESS MATRIX FOR THE GLOBAL DISPLACEMENT 
SOLUTION. INCREASE FIELD LENGTH BY (OCTAL) 

THE STIFFNESS MATRIX CANNOT BE DECOMPOSED 


INSUFFICIENT BLANK COMMON AREA IS AVAILABLE TO STORE THE LOAD 
VECTOR FOR THERMAL STRESS CALCULATIONS. INCREASE FIELD LENGTH 
BY (OCTAL). 
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LCNTRL 


Double root could not be successfully bypassed. 


901 

902 Zero determinant could not be successfully bypassed. 

903 Search for buckling condition did not stabilize in specified 
number of iterations. 

Check for panel being in buckled condition for zero load. 

UPRBND 

BUCKLED CONDITION element ID, subelement, value (calculated) 
value (input) 

Input load on subelement caused buckled condition 
62 LOAD FIELD IS ALL ZERO OR ALL TENSION 

64 PANEL WAS DETECTED TO BE IN A BUCKLED CONDITION WITH REPECT 

TO THE INITIAL LOAD DISTRIBUTION 
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DM 

30 

GAPUPP 

DB 

906 

SBLKDT 


DBLERT 

STORE 

22 

24 

PRDMTX 

5 

PROOT 


Double root encountered in upper bound calculation 
Matrix reduction failed 
Matrix block definition faulty 
The MATRIX HAS A ZERO ROW 

MATRIX SIZES OR RELATIVE POSITIONS ARE INCOMPATIBLE 
Zero Determinant at Block xxx 

Double Root Detected 

Node number not found in substructure row definitions 
Node number not found in substructure column definitions 


Element nodal interconnective matrix singular 


ZARK failed to converge in the maximum number of iterations 
specified 
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PROOT 

ZARK failed - A zero in the path of a subsequent one. 
Complex root found that is not one of a conjugate pair 

SOLVV 

8 SINGULARITY ENCOUNTERED IN ATTEMPT TO PRODUCE STIFFNESS 

MATRIX 

DISPLAC 

INSUFFICIENT CORE IS AVAILABLE TO CONSTRUCT THE STIFFNESS 
MATRIX FOR THE EIGENVECTOR SOLUTION INCREASE FIELD LENGTH 
BY 

THE STIFFNESS MATRIX CANNOT BE DECOMPOSED 

FETCHD 

14 Attempt fetch a non-existent nodal displacement component 

READTR 

10 ID on nodal transformation invalid 

12 ID on elemental transformation invalid 

FNLDIS 

20 Attempt to invert elemental tranformati on matrix failed. 
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5.0 RESTRICTIONS 

The following restrictions apply to this version of the program. 

5.1 Analysis Oriented Restrictions 

a. Prebuckling deformations are ignored and linear theory (small 
deformations and linear elastic materials) is used. 

b. Each layer is orthotropic with respect to plate axes. 

c. Plate elements have constant thicknesses. 

d. The cross-section is uniform in the axial direction. 

e. Loaded edges of each element forming the cross-section are simply 
supported. 

f. The arc of a curved plate element has to be less than or equal to 
180 degrees. 

g. Buckling load will not be computed for a structure already buckled 
by the input loads. 

h. The temperature distribution in the longitudinal direction is a 
constant. 

i. Only plate in-plane biaxial loads and beam axial loads are considered 
as applied loads in the buckling formulation. Further, they are 
treated as constant values. 

j. For the input loads option, plate biaxial loads are constants in 
each subelement of a plate element. 
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5.2 Programming Oriented Restrictions 


a. Maximum number of elements definable in three substructures is 35, 
and the maximum number of nodes is 35. 

b. Maximum number of types of plate elements is 10. These are of the 
same type in the sense that their lamina stiffness matrices, are 
of the same content. 

c. The maximum number of layers is 25 for plate elements. 

d. Maximum number of beam elements is 30. 

e. The maximum number of layers is 35 for beam elements. 

f. Two real roots of the determinant expression det (DT) = 0 are 
considered double if they differ by less than 0.00003%. 

g. The imaginary part of the complex roots of the determinant 

expression det (DT) = 0 is exactly set to zero if its numerical 

-6 -5 

value is less than 10 , or when it is less than 10 times the 

real part of the number. A similar test applies to the real part 

of the number. 

h. It is assumed that no coupling exists between bending and stretching 
when all of the elements of the B-matrix are less than 1.0. 

i. Maximum number of wave numbers that can be investigated in one data 
set is 30. 

j. The minimum number of elements is either 2 plates or 1 plate and 
1 beam. 

k. There must be at least 1 interior node in each substructure type. 

l. Node restrictions in substructues ; see Section 2.5. 

m. An element cannot have nodes both in the start and end substructures, 
which means closed cylinder problems must be done with a single 
substructure. 
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n. Number of buckling subdivisions per plate is between 2 and 20. 

o. Total number of elements definable with thermal or external 
loading is 100. 

p. Each substructure is restricted to 14 elements. 

q. Thermal stress calculation is limited to twelve points in the 
x-di recti on and is set to five in the y-di recti on. 
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5.3 Mathematical Oriented Restrictions 

In the buckling calculation the process is done strictly incore where 
the number of repeat substructures is arbitrary, but will be limited 
in practice by the accumulation of rounding errors. 

The eigenvector computation requires that the buckling load be well 
isolated. Problems which buckle without a determinant change of 
sign will have multiple eigenvectors for the buckling load. 

The program appears to have some modelling restrictions relative 
to use of high wave numbers in buckling calculations. See appendix 
for details. 
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6.0 PROGRAM INPUT 


The data input to this program consists only of cards, and no data tapes 
are required. 

6.1 Input Data Format 

Fields of 5 or less are right adjusted integers. Fields of 10 are floating 
point numbers with format F10.0. 

CARD 1 (8A10) 

Columns Column 1 must be blank. 

2-80 Title of run. Any characters anywhere in columns. This title 

is printed out in several strategic places in the output for 
the purpose of identification. 

CARD 2 (1615) Intermediate Print Control (IPC) Array 

The intent of the options on this card is for diagnostic checking. 

Col umns 

I- 5 I PC ( 1 ) = 1 Print intermediate results (1) 

= Blank Suppress print 

This option prints out control lists obtained from input data. 

6-10 I PC ( 2 ) = 1 Print intermediate results (2) 

= Blank Suppress print 

This option prints out the p roots of the determinant 
expression det (DT) = 0 (equilibrium equations) during the 
upper bound calculation. 

II- 15 I PC ( 3 ) = 1 Print intermediate results (3) 

= 31ank Suppress print 

This option prints out the elemental and merged stiffness 
matrices along with the determinant calculation information 
during the upper bound calculation. 
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Columns 


16-20 IPC (4) = 1 Print intermediate results (4) 

= Blank Suppress print 

This option prints out information about the element upper 
bound initial estimates and the actual upper bound. 

21-25 I PD ( 5 ) = 1 Print intermediate results (5) 

= Blank Suppress print 

This option prints out the p roots of the determinant 
expression det (DT) = 0 (equilibrium equations) during the 
panel buckling calculation. 

26-30 I PC (6) = 1 Print intermediate results (6) 

= Blank Suppress print 

This option prints out the elemental matrices during the 
panel buckling calculation. 

31-35 I PC ( 7 ) = 1 Print intermediate results (7) 

= Blank Suppress print 

This option prints out the merged stiffness matrices and 
determinant calculation information during the panel buckling 
calculation. 

36-40 IPC (8) = 1 Print intermediate results (8) 

= Blank Suppress print 

This option prints out the compacted stiffness matrix, the 
element p roots and the element to. . 

41-45 IPC ( 9 ) = 1 Print intermediate results (9) 

= Blank Suppress print 

This option prints out the p roots of the determinant 
expression det (DT) = 0 (equilibrium equations) during the 
thermal stress calculation. 
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Columns 


46-50 I PC (10) = 1 Print intermediate results (10) 

= Blank Suppress print 

This option prints out the element matrices during the 
thermal stress phase. 

51-55 I PC ( 1 1 ) = 1 Print intermediate results (11) 

= Blank Suppress print 

This option prints out the merged stiffness matrix. 

56-60 IPC(12) = 1 Print intermediate results (12) 

= Blank Suppress print 

This option prints the merged right-hand side vector, the 
calculated global deflection vector and other check prints 
in the deflection calculation. 

61-65 I PC ( 1 3 ) = 1 Print intermediate results (13) 

= Blank Suppress print 

This option turns on the check prints in the stress 
calculation. 

66-70 I PC (14) = 1 Print intermediate results (14) 

= Blank Suppress print 

This option prints out subelement force and deflection 
matrices during buckling calculation 
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CARD 3 (1615) Program Control Array 

Col umns 

I- 5 JPC(l) = 1 Calculate panel upper bound only 

= 0 (Blank) Proceed as normal 

6-10 J PC ( 2 ) = 0 (Blank) No relative displacements 

= 1 Compute relative displacements 

II- 15 JPC(3) = 1 Data check only is performed 

= 0 (Blank) Non data check only mode 

16-20 JPC(4) n Number of iterations allowed in buckling 

calculation, where 1 < n < 100 
= Blank Number will be defaulted to 100 

21-25 JPC(5) n Number of iterations used in eigenvector 

calculation, where 4 <_ n ^ 100 
= Blank Number will be defaulted to 4 

Option should be used with care. Can be used to find 
multiple eigenvectors for problems with coincident roots by 
running problem with n iterations and then with n + 1 
iterations. 

26-30 JPC(6) Default plate element thermal input 

option 

= 1 Input Tp Z (constant) at mid-plane of 

each lamina 

= 2 Input element coefficients A z , B z> C z 

2 

of the equation T pz = A z + B z z + C_,z 

for each element, z is along the local 
2 axis of the element and its origin is 
at the outer surface of the element. 
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Columns 


31-35 


36-40 


41-45 


46-50 


JPC( 7) 


JPC(8) 

JPC(9) 


JPC(10) 


= NUMX 


- NUMX 


Default value for number of plate element 
subelements used in buckling analysis. 
Value must be greater than 1 and less 
than or equal 20. To specify number of 
subelements other than default value, see 
Card P. 

= First Fourier harmonic (odd number) for 
output of stresses and displacements 

= Last Fourier harmonic (odd number) for 
output of stresses and displacements. 

(The stresses and displacements are 
summed from 1 and are printed for 
harmonics between specified limits) 

Output option for stresses and displace- 
ments 

where NUMX = number of equally spaced 
points along the x-axis for the output of 
stresses and displacements 
Table of NUMX values for output of 
stresses and displacements will be given 
on Card 8 (To be used when output at 
points unequally spaced along the x-axis 
is desired) 



CARD 4 


Problem Characteristics 


Col umns 



1-5 

Number of Nodes in the section (i.e., in start substructure, one 
repetitive substructure and the end substructure) 

6-10 

Number of Beam elements in the section 


11-15 

Number of Flat Plate elements in the section 


16-20 

Number of Curved Plate elements in the section 


41-50 

Length of the section 


51-60 

Transverse Load or Strain value 


61-70 

Biaxial Load Ratio for first plate element 
Input only when this option is selected. 


71-75 

Maximum Fourier harmonic used to describe constant thermal 
along local x-axis (Use only odd numbers). 

load 

For the transverse load or strain value, input a negative value for tension 
and a positive value for compression. Biaxial load ratios are computed with 
the same sign convention as above. 
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CARD 5 
Columns 

I- 5 

= 4 
= 5 

= 6 

= 7 

= 8 

6-10 MOPT 
= 1 

= 2 

v = 3 

= 4 

II- 15 MM I 

16-20 MMA 


Analysis Options 


Load/Stress Options 

Calculate Thermal Stresses only (see Section 1.1). 

Buckling due to axial load under a thermal environment 
(see Section 1.2). 

Buckling due to thermal effects only. Critical temperature 
ratio is computed (see Section 1.2). 

Buckling due to (force/length) for a panel that is 
initially prestressed by specified biaxial loads. (See 
Section 1.3). 

Buckling analysis, where the critical load ratio is 
computed. See Section 1.3. 

Option control for the loop on the longitudinal wave number M. 
Start the loop at wave number MMI and loop until a minimum 
load is found, then interrupt (Max. 30 loops) 

Start the loop at wave number MMI and loop to wave number 
MMA (Max. 30 loops) 

Start the loop at first value of M-list and loop through 
M- list until a minimum load is found, then interrupt (M-list 
must be input on Cards 6,7) (Max. 30 loops) 

Start the "bop at first value of M-list and loop through 
M-list. (M-list must be input on Cards 6,7) Max. 30 loops) 

Starting value for the loop on the longitudinal buckling wave 
number M 

Set = 1 for options 3 and 4 

End value for the loop on the longitudinal buckling wave 
number M 

Set = number of values in M-list for options 3 and 4 
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Columns 

21-25 NLW Lower Limit for root searching criteria 
(If left blank, NLW is defaulted to 0) 

26-30 NUP Upper Limit for root searching criteria. If left blank, 

NUP is defaulted to 1 NLW >_0, NLW < NUP. The program does 
a complete solution of the data for each root searching 
criteria n/n+1 where NLW <_ n < NUP. (n is the number of 
roots below the trial load.) 

31-35 Default number of element subdivisions for relative displace- 

ment calculation (defaulted to 5 , if left blank) 

41-50 SLW Lower bound (first plate element) for search interval (defaulted 
to 0 . 0 , if left blank) 

51-60 SUP Upper bound (first plate element) for search interval. If left 
blank an upper bound is calculated for the structure for each M 
value. If specified SLW SUP. Note that the same SLW, SUP 
pair is used for every M value to be searched. 

See Section 1 for discussion of load options 4 to 8 . Refer to Section 2.6 

for discussion of upper bounds. For a discussion of the root searching 

criteria, see Section 6.1 of Reference 2 0 
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CARD 6 

M- list (limited to 30 values) 
(List of wave numbers) 

Columns 


1-5 

First value of M-list 

6-10 

• 

» 

Second value of M-list 

• 

• 

75-80 

Sixteenth value of M-list 

CARD 7 

M-list (continued) 

(Use this card if necessary; if not omit) 

1-5 

• 

• 

Seventeenth value of M-list 

• 

• 

65-70 

Thirtieth value of M-list 
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CARD 8 

Columns 

1-10 

11-20 

m 

71-80 

CARD 9 
Columns 
1-10 
11-20 
21-30 
31-40 


Table input for output of stresses and displacements 
Omit if JPC(IO) >_ 0 or blank 


First x-coordinate 
Second x-coordinate 


Eigth x-coordinate 


Ninth x-coordinate 
Tenth x-coordinate 
Eleventh x-coordinate 
Twelfth x-coordinate 
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CARD T-l 


Lamina Thickness Table 


Columns 

I- 5 Letters "THICK" 

II- 80 , value, .../ 

where value may be + nnn.nnn 
or + n.nnnn E nn 

Number of significant digits is limited to 14 

The table is terminated by the character slash (/). If the entries of a 
table exceed one card, the table may be put on more cards using columns 
11-80 (columns 1-10 blank) with the last table value followed by a slash 
The values in the table are delimited by commas and physical end of card 

Omit card if Table option not used for element property input. See 
Section 6.4 for discussion of table use. 

Note: Only one thickness table is allowed per data set. 
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CARD T-2 

Lamina Material Properties Table 

Col umns 


1-5 

Letters "TABLE" 

9-10 

Table Number 

11-80 

...» ... ) value, . . . / 


Same rules for this table as the one on CARD T-l 

Omit card if Table option is not used for element property input. 
See Section 6.4 for discussion of table use. 


Table Lengths 
El ement Length 

Plate 4 

Rectangular 

Beam 3 

Ci rcular 

Beam 2 
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CARD C Substructure Interrelationship Element Pair 

Columns 


1-1 


3-5 

8-10 


73-75 


Letter "C" 
Element i 


Pair 1 (this card) 


Element j J 

Element i of start-repeat substructure is interrelated to 
Element j of repeat-end substructure 

Element i 1 


Pair 8 (this card) 


78-80 Element j J 

Repeat this card for each set of 8 interrelationship element pairs until all 
interrelationship element pairs have been input. 

Note: Interrelationship elements must be included if the number of repeat 

blocks > 2. See Section 2o5.2 for discussion. 
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CARD E 


Equivalent Nodes 


Columns 

1-1 

3-5 

8-10 


Letter "E" 


Node i 
Node j 


Pair 1 (this card) 


73-75 

79-80 


Node i 
Node j 


Pair 8 (this card) 


Repeat this card for each set of 8 equivalent node pairs until all 
equivalent node pairs have been input. See section 6.2 for 
discussion of equivalent nodes. 

Note : Equivalent nodes must be interior nodes. This card is used to 

specify those nodes which move together structurally. Omit this 
card if no equivalent nodes exist. 
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CARD N 
Columns 


1-1 


3-5 

10-10 


21-30 

31-40 


Nodal Definition 


Letter "N" 

User Node Number 
Nodal Boundary Condition 

= Blank Interior 

1 Simple Support 

2 Clamped 

3 Free 

4 Sprung (Card S must follow) 

Y - coordinate of node in right-handed global coordinate 
system 

Z - coordinate of node in right-handed global coordinate 
system 

Internal node number is determined by input position, i.e., 
first node input is internal node one, etc. 
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CARD D 


Panel Substructure Definition 


Columns 

1-1 Letter "D" 

3-5 The number of substructures that the buckling determinant is 

divided into. This includes the start and end substructures 
plus the repeat substructures. 

8-10 User node number which the Lowest Internal Node Number in 

Start Substructure 

13-15 User node number which is the Highest Internal Node Number 

in Start Substructure 

18-20 User node number which is the Lowest Internal Node Number in 

Repeat Substructure 

23-25 User node number which is the Highest Internal Node Number 

in Repeat Substructure 

28-30 User node number which is the Lowest Internal Node Number in 

End Substructure 

33-35 User node number which is the Highest Internal Node Number 

in End Substructure 

A substructure is left undefined by not specifying substructure definition 
node pair. 

Permissible Combi nations 

Start Substructure Only 

Start and End Substructures Only 

Start, Repeat (one or more) and End Substructures 

It is recommended that the start substructure be larger of the 

three substructures, for efficiency reasons in the determinant 

calculations. Note: See card N for definition of Internal Node 

Number. 
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CARD S Nodal Spring Stiffness 

Columns 

I- 1 Letter "S" 

3-5 User Node Number 

II- 20 Local w component for node 

21-30 Local 9 component for node 

31-40 Local v component for node 

41-50 Local u component for node 

See Section 2.4 for discussion of sprung nodes. 
Omit this card if there are no sprung nodes. 
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CARD F 
Columns 
1-1 
3-5 

11-20 0FF1 


21-30 0FF2 

31-40 0FF3 
41-50 0FF4 


Plate Element Offset 


Letter "F" 

User Plate Number 

Offset z in the local z direction at the starting end 
(y = -b/2) of the plate element, measured positive in the 
positive z direction, from the negative surface of the 
element to the grid. 

Offset y Q in the local y direction at the starting end 
(y = -b/2) of the element, measured positive in the positive 
y direction, from the end of the element to the grid. 

Offset z Q at the end y = +b/2 of the plate element. 

Measured similarly to 0FF1. 

Offset y Q at the end y = +b/2 of the plate element. 

Measured similarly to 0FF2. 

Omit this if no plate offsets exist. 

See Section 2.3 for discussion of offsets. 
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CARD P 
Columns 
1-1 
3-5 
8-10 
13-15 
20-20 


25-25 


29-30 

35-35 


40-40 


43-45 


50-50 


Plate Definition 


Letter "P" 

User Plate Number 

Plate definition Node I 

Plate definition Node J 

Plate element type 
= 1 Flat Plate 
2 Curved Plate 

Element Offset Switch 
= 0 (Blank) No offsets for this plate 
1 Element has offsets 

Number of laminas in element 

Input type for element material properties 

= 0 (Blank) P-1 or P-2 cards are used for material property 
input 

1 Material properties are input through tables 

N ^2 1 oad selection switch 

= 0 (Blank) The N^ load is applied this plate ("affected") 
1 This plate will have an effective l° ac * ^ 

ZERO (Available for options 1,2,3) 

Number of subdivisions used in relative displacement calcu- 
lation. If left blank, the default number specified on 
Card 5 is used. 

Input option for material properties 

= 0 (Blank) Engineering constants E-j-j, E^* ^12 are 

specified 

1 Lamina stress-strain matrix Q is specified 
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Columns 


51-55 


56-60 


61-65 


66-70 

71-80 


Local value for plate thermal input (Required for LOADOP = 4,5,6) 
= 0 Default option is used 

1 Input T z at mid-plane of each lamina (constant) 

2 Input element coefficients A z , B z> C z 

3 No thermal input 

Load distribution input option (Required for LOADOP = 7,8) 

= 0 Load input present 
1 No load input 

Substructure to which the element is assigned. If both nodes 
are in a substructure, then the element is in that substructure, 
otherwise: 

If one node is in the start substructure and one node is in 
the repeat substructure, then the element is in the repeat 
substructure. 

If one node is in the repeat substructure and one node is in 
the end substructure, then the element is in the end sub- 
structure. 

If one node is in the start substructure and one node is in 
the end substructure, then the element is in the end sub- 
structure. 

(Required for LOADOP =4-8) 

1 = Start substructure 

2 = Repeat substructure 

3 = End substructure 

Number of subelements used in buckling calculation. (If 
left blank, default value is used.) 

Curved Plate Element Radius (see Section 2.2.1) 
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CARD P-1 


Plate Thickness and Material Properties 


Columhs 



1-10 

T 

Thickness of lamina 

11-20 

E 11 

E- modulus for direction 1 

21-30 

E 22 

E- modulus for direction 2 
isotropic laminas. 

31-40 

RNUA 

Poisson's ratio, 

41-50 

G 12 

G- modulus 


E ^2 neec * n °t be entered for 


The subscript 1 denotes the longitudinal axis and the subscript 
2 the transverse axis of the plate local coordinate system. 
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CARD P-2 
Col umns 

Plate Lamina Thickness 

and Q. . Matrix 
U 

1-10 

Thickness of lamina 


11-20 

Qll element of lamina 

Q. . matrix 
1 J 

21-30 

Q ^2 element of lamina 

Q. . matrix 

* J 

31-40 

Q ^2 element of lamina 

Q.. matrix 
1 J 

41-50 

Qgg element of lamina 

Q. . matrix 

' J 


Q. . are the elements 
U 

of the stress-strain matrix 


equation (4.3) of Reference 2. 
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CARD P-3 Plate Temperature Distribution Coefficients 

Columns 

I- 10 TX Tp x Constant value to describe the temperature 

distribution in the x-direction. 

II- 20 TYA a y for T py = a y + b y .y + c y .y 2 , 

21-30 TYB b y ► temperature variation in the local y-directiort 
31-40 TYC c 

7 ^ - - 2 
41-50 TZA A z for T pz = A z + B z z + C.,z » 

51-60 TZB B z , temperature variation in the local z-di recti on. 

61-70 TZC C z 

J 

(Omit columns 41 to 70 if T pz is input at mid-plane) 

The P-3 card follows the P-1 or P-2 (or equivalent) for a plate element. 
For elements in the repeat substructure, the sequence of P-3 and P-4 cards 
must be repeated for each repeat substructure. 
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CARD P-4 


Plate Lamina Thermal Properties 


Columns 

I- 10 ALPHX 

II- 20 ALPHY 

21-30 TZ 

(Omit columns 
There must be 


Linear Thermal Coefficient of Expansion in the 
x direction 

Linear Thermal Coefficient of Expansion in the local 
y direction 

Value for T 

21-30 if A , B , C z values are input) 
a P-4 card input for each lamina in plate. 
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CARD P-5 N^-j Load Distribution Input 
Columns 

I- 10 N-J.J value for first subelement 

II- 20 N-j i value for second subelement 

21-30 N.ji value for third subelement 

31-40 N-ji value for fourth subelement 

41-50 N^i value for fifth subelement 

51-60 N-ji value for sixth subelement 

61-70 N-j i value for seventh subelement 

71-80 N-ji value for eigth subelement 

For elements with more than eight subelements, repeat the P-5 card until 
all N -j -j values are input for all subelements. (Limit of 20 subelements 
per element) 

For elements in the repeat substructure, there must be a group of P-5 cards 
for each repeat substructure. 

The tensile loads are input as positive values and the compressive loads 
as negative values. 
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CARD P-6 N^2 Load Distribution Input 

Columns 

I- 10 N22 value for first subelement 

II- 20 ^2 value for second subelement 

21-30 N22 value for third subelement 

31-40 N22 value for fourth subelement 

41-50 N22 value for fifth subelement 

51-60 N22 value for sixth subelement 

61-70 N22 value for seventh subelement 

71-80 N22 value for eigth subelement 

For elements with more than eight subelements, repeat the P-6 card until all 
N -j 1 values are input for all subelements. (Limit of 20 subelements per 
element) 

For elements in the repeat substructure, there must be a group for P-6 cards 
for each repeat substructure. 

Tensile loads are input as positive values; compressive loads as negative 
values. 
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CARD B 
Columns 


1-1 

3-5 

8-10 


13-15 


18-20 

25-25 


26-30 


31-35 


Beam Definition 


Letter "B" 

User Beam Number 
Beam Definition Node 
Type of Beam 
= 1 General 

2 Rectangular (may be laminated) 

3 Circular (may be laminated) 

Number of layers if laminated 

Input type for element material properties 

= 0 (Blank) B-l , B-2 or B-3 cards are used for material 
property input 

1 Material properties are input through tables 

Thermal Input Option 

= 0 Thermal Data is input for this element (i.e., a 

B-4, B-5, B-6, or B-7 card is input) 

1 No Thermal Data 

Load Input Option 

= 0 Loads are input for this element (i.e., a B-8 

card is input) 

1 No load data 


36-40 Substructure to which the element is assigned (Required 

for L0AD0P = 4-8. See P card for explanation.) 
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41-50 

51-60 


Angle between the local y-axis of beam element and the global 
y-axis (measured clockwise, from the global y-axis) 

Area of beam (Option 1) 

Beam element definition node must be equivalenced to a node 
of plate element. 
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CARD B-l 


General Beam Properties 


Columns 

I- 10 E-modulus of beam element in longitudinal direction 

II- 20 G-modulus of beam element material 

21-30 Moment of inertia about local y-axis 

31-40 Moment of inertia about local z-axis 

41-50 Warping constant for beam element about shear center 

5 1 - GO Torsion constant of bead or lip about shear center 

61-70 y distance measured from the shear center to the centroid of the 

beam (measured parallel to local y axis) 

71-80 z distance measured from the shear center to the centroid of the 

beam (measured parallel to local z axis) 

See Section 2.2.2 for description of local coordinate system. 
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CARD B-2 


Rectangular Beam Lamina Thickness and Material Properties 


Columns 

I- 10 Lamina Thickness 

II- 20 E-modulus of lamina in longitudinal direction 

21-30 G-modulus of lamina material 

31-40 Width of beam element 

If the width of all lamina is the same, then only the width 
for the first lamina need be specified. 
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CARD B-3 


Circular Beam Lamina Radius, and Material Properties 


Col umns 

I- 10 Lamina Outer Radius 

II- 20 E-modulus of lamina in longitudinal direction 

21-30 G-modulus of lamina material 

For a tube element, input the hole as the first lamina with zero E and G. 


m 
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CARD B-4 


General Beam Thermal Loads 


Columns 


1-10 
11-20 
21-30 

For elements in the repeat substructure, there must be a B-4 
for each repeat substructure for the element. 


MTY 

Thermal 

Moment 

about 

the 

local 

y axis, M yT 

MTZ 

Thermal 

Moment 

about 

the 

local 

z axis, M^ 

PIT 

Thermal 

axial 

load - 

P T 




The B-4 card follows the B-l card for a beam. 
See equation (4.53) to (4.55) of Reference 2. 


card input 


}S , 


6.32 



CARD B-5 


Rectangular Beam Temperature Distribution Coefficients 


Col umns 


1-10 

TXO 

T bx 

constant value to describe the temperature 
distribution in the x-direction 

11-20 

TAZ 

a 

z 


21-30 

TBZ 

b 

z 

for T, = a + b'z + c *z , 

, bz z z z 

temperature variation in the local 

31-40 

TCZ 

c z 

J 

z-di recti on 

The B-5 card 

follows the B-2 (or equivalent) for a beam. For elements in 


the repeat substructures, the sequence of B-5 and B-6 cards must be repeated 
for each repeat substructure. 
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CARD B-6 


Rectangular Beam Lamina Thermal Properties 


Columns 

I- 10 ALPHX Linear Thermal Coefficient of Expansion in the x-direction 

II- 20 TY T, - constant value to describe temperature distribution 

by 

in the local y-di recti on 

There must be a B-6 card input for each lamina of the beam. The B-6 card(s) 
for a beam follow the beam's B-5 card. 
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CARD B-7 


Circular Beam Lamina Thermal Properties 


Col umns 

I- 10 ALPHX Linear Thermal Coefficient of Expansion in the x-di recti on 

II- 20 TR Lamina Temperature 

The B-7 card(s) for a beam for the B-3 (or equivalent) for the beam. There 
must be a B-7 card input for each lamina of the beam. 

For elements in the repeat substructures , there must be input a group of 
B-7 cards for each repeat substructure. 
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CARD B-8 


Beam Load Input 


Col umns 

1-10 Axial load for beam, tension = (+) and compression = (-). 

The B-8 card follows the B-l, B-2 (or equivalent), or B-3 (or equivalent) 
for a beam. For elements in the repeat substructure, there must be a B-8 
card input for each repeat substructure. 


x 
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CARD L Last Card of Data Set 

Col umns 

1-1 Letter "L" 

This card must terminate each data set. 
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6.2 Use of Equivalent Nodes 

Equivalent nodes are used to cause two nodes to move together 
structurally. Beam are attached to plate elements so beam definition 
nodes must always be made equivalent to plate nodes. Care must be 
exercised so that (1) the use of equivalent nodes does not 
alter a substructure definition, and (2) a node that is a left 
(right) node of an equivalent pair does not appear as a right 
(left) node of a subsequent equivalent pair. See examples below: 

not permissable permissable 
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6.3 Data Stacking 

LOADOP = 1, 2, 3 


ORDER 

1 

2 

3 

4 

5 

Next if required 


CARD 
1 
2 

3 

4 

5 

6 
7 
N 
D 

T-l 
T-2 

C 
E 
S 
F 
P 

P-1 (or equivalent) 
or (not both) 

P-2 (or equivalent) 
B 

B-l 

B-2 (or equivalent) 
B-3 (or equivalent) 
L 


Omit any cards that are inapplicable 


Repeat for 
each lamina 


Repeat for 
each lamina 


All plate definition cards 


All beam definition cards 
(Omit if no beams) 

Last card in data set 


Repeat this order of card stacking for multiple data sets. 
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LOADOP = 4 


CARD 

1 

2 

3 

4 

5 
8 
9 
N 


ORDER 

1 

2 

3 

4 

5 

Next if required 


or 


or 


D 

T-l 




T-2 


C 

E 

S 

F 

P 

P-1 

P-2 

P-3 

P-4 

B 


I Omit any cards that are inapplicable. 


(or equivalent) 
(or equivalent) 

Repeat for each 
lamina 



Repeat for each lamina 

Repeat for elements in 
Repeat Substructure 


B-l 


) 


j 


PI ate 

Definition 

Cards 




B-4 


B-2 (or equivalent) Repeat for each lamina 

B-5 I Repeat for elements 

B-6 Repeat for each lamina | in Repeat Substructure 

B-3 (or equivalent) Repeat for each lamina 

B-7 Repeat for each lamina I Repeat for elements in 

J Repeat Substructure 


All beam 
defini tion 
Cards (Omit 
if no beams) 
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LOADOP = 5 or 6 


CARD 

1 

2 

3 

4 

5 

6 
7 
N 


ORDER 

1 

2 

3 

4 

5 

Next if required 


D 

T-l 

T-2 


t 

E 
S 
F 
P 

P-1 (or equivalent) 
P-2 (or equivalent) 
P-3 


Omit any cards that are inapplicable 


or Repeat for each lamina 

( Repeat for elements 
in Repeat Substructure 


I 


Plate Definition 
Chart 


B 


B-l 

B-4 


or f 


B-2 (or equivalent) 
B-5 


B-6 Repeat for each 
B-3 (or equivalent) 
B-7 Repeat for each 


Repeat for each lamina 

1 Repeat for elements in 
Repeat Substructure 

r > 

Repeat for each lamina 

lamina 1 Repeat for elements 

J in Repeat Substructure 


A1 1 beam definition 
Cards (Omit if no 
beams) 
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LOADOP = 7 or 8 


CARD 

1 

2 

3 

4 

5 

6 
7 
N 
D 

T-l 

1-2 

C 

E 

S 

F 

P 

P-1 

P-2 

P-5 

P-6 

B 

B-l 


ORDER 

1 

2 

3 

4 

5 

Next if required 


Omit any cards that are inapplicable 


or 


or 


B-2 

B-3 

B-8 


(or equivalent) I Qr Repeat for Each Lamina 
(or equivalent) J 

( Repeat for elements 
in Repeat Substructure 


(or equivalent) 

(or equivalent) 

Repeat for elements in Repeat Substructure 


Plate Definition 
Cards 

All beam definition 
cards (Omit if no 
beams) 
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6.4 Use of Tables 


The program has a table input option to facilitate the repetitive 
nature of the input for the thicknesses, radii and material properties 
of the plate and beam elements. This option replaces the P-1, P-2 
B-2 or B-3 cards. The first item of the P-1, P-2, B-2 and B-3 cards 
(lamina thickness or lamina radius) is specified on the T-l card. 

All other items of the P-1, P-2, B-2 and B-3 cards are input on the 
T-2 card. These items are the material properties, E, v, G or Q, 
for each layer. Each set of values on the T-2 card will be identified 
by a user table number. In a data set, it is permissible to use both 
input options for different elements, but the same input option must 
be used for all layers of any one element. 

To use this table input, first input the proper parameter on columns 
35 and 25 of the P and B cards, respectively. Then, on subsequent 
cards, one or more for each element, specify the table numbers for 
element material properties and the locations of the thicknesses or 
radii on the T-l card. An example of the card is shown below. 

(1,2,3,,,,!/ 3,2, ,5, 6, 7, 6) 

thickness material 

or radii properties 

The left and right parenthesis are used to initiate and terminate the 
data for one element. The / is used to separate the T-l and T-2 
entries. The numbers to the left of / corresponds to the n*"* 1 entries 
of the T-l card. The numbers are input in the order of the layers 
and separated with commas. If a number is omitted between commas, 
the previous integer is used for that layer. The numbers of commas 
to the left and right of / must always be equal. 
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The data for the example above will be interpreted in the following 
manner. 


Layer No , 


n** 1 Item of T-1 Card 


Table No. of T-2 


1 

2 

3 

4 

5 

6 
7 


1 

2 

3 

3 

3 

3 

1 


3 

2 

2 

5 

6 
7 
6 
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7.0 OUTPUT 


The printing of the input data is coordinated with the description of 
the input in Section 6.1. The items of output that are produced are 
as follows: 

1 . Header Page 

2. Program Control Options 
Card 3 

3. Problem Characteristics 
Card 4 

4. Analysis Options 
Cards 5, 6 and 7 

5. Nodal Data 
Card N 

6. Substructure Definition Data 
Card D 

7. Repeat Substructure Interrelationship Data 

Card C. Element I connects the start-repeat substructure and 
Element J connects the repeat-end substructure. 

8. Equivalent Node Data 

Card E Node I is equivalent to Node J. 

9. Table Data 
Cards T1 and T2 



10. Element Data 
Plates 

Cards P, P-1 (or equivalent) or P-2 (or equivalent) 

Each plate is accompanied by its A-, B- and D-matrices. The 

A-matrix represents the extensional stiffness of the plate 

element, while the B- and D-matrices are the coupling stiffness 

and bending stiffness, respectively. For thermal problems cards 

P-3 and P-4 will be printed and for input loads problems, P-5 and P-6. 

The A, B and D matrices are calculated relative to the element midplane. 

Beams 

Cards B, B-l, B-2 (or equivalent) or B-3 (or equivalent) 

For thermal problems, either cards B-4 or B-5 and B-6 or B-7 will 
be printed. For input loads problems, B-8 card will be printed. 

11. Plate Element Input Offsets 
Card F 

12. Plate Element Net Offsets 

The net Zo offsets are measured from the midplane of the element 
to the point that is offset to. 

13. Element Transformation Data 

The plate element widths, the sines and cosines of the transformation 
angles (local to global) at the definition nodes and the element 
types are printed. The sign convention for the angles are shown 
in Figure 5.2 of Reference 2. The plate width for curved elements 
is the arc length between nodes. 

14. Element Substructure Assignment Data 
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15. Field Length Data 

The field length for the buckling calculation is exact, whereas 
the one for the eigenvector and thermal stresses calculation is 
an upper bound to the one needed. 

16. Thermal Stresses and Displacements 

The element thermal stresses and displacements are output with 
respect to the local coordinate systems. For every specified 
x-value, these stresses and displacements are printed at five 
equidistant local y-values of the plate element (LOADOP = 4). For 
LOADOP = 5, the and average loads are printed for each 
subelement when requested. (This is accomplished by specifying 
first and last print Harmonic. See Card 3, fields 8 and 9.) 

17. Panel Upper Bound 

18. Buckling Load or Critical Ratio Search History 

The number of negative elements found on the diagonal of the 
panel matrix during the determinant calculation. The buckling 

D 

determinant = A.(2 ). There is a history printed for each wave 
number. This number of negative elements indicates the number of 
eigenvalues that is below the trial load. For LOADOP = 6 and 8 
(critical ratio options), the ratios are printed instead of trial 
loads. 

19. Element Loads 

The inplane biaxial loads, and an ^ the biaxial strains, 

E-|i and E^£ f° r plate elements are printed here. For beam elements, 
the axial loads and strains are printed. The total loads, axial 
and transverse are also included. Sign convention used here is 
such that tension is (-) and compression is (+). For LOADOP = 6 or 8, 
there will be no element loads print. 
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20. Buckling Load Computation Summary 

The critical wave number among those searched is identified. The 
buckling load for the first element of the panel is printed for 
each wave number. For L0AD0P = 6 and 8, the critical ratios are 
summarized. 

21. Eigenvector 

For each iteration, the eigenvectors for every freedom of all nodes, 
with the exception of constrained or sprung nodes, are printed along 
with the normalizing factors. User should inspect these eigenvectors 
to check on convergence. These values are printed with respect to 
the global coordinates. 

22. Relative Displacements 

The element displacements are relative to the element local 
coordinate system. 
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8.0 SAMPLE PROBLEMS 


The input, output and accompanying sketches of four (4) sample 
problems are included in this section. The circled items on the 
figures and discussions are user element numbers and the uncircled 
integers are user node numbers „ 

8.1 Example 1: Potpourri Input Problem 

The sample input for this problem is presented for this rather 
arbitrary structure, shown in Figure 8.1. This example is included 
here to illustrate many input features of BUCLASP3. 

The total structure is assumed to comprise of two more repeat 
substructures than is shown in Figure 8.1. .Node 101 is simply 
supported, node 117 clamped, node 107 sprung, and nodes 104 and 
113 are free to deform. Plates (f^) and (f§) are separated by an 
offset of 1 inch. ([§) is a three layered circular beam and (8) 
is a general beam element that is offset 2 inches from the midplanes 
of plate elements (6) and ( 9 ) . 

Each of the plate elements (T) through ( 7 ) is subdivided into 
3 subelements. Plate element (f8) is subdivided into 5 subelements 
and each of the remaining elements is subdivided into 2 subelements. 
The critical load ratio option is used here and the plate subelement 
and beam input loads are specified. 
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Repeat 

Substructure 


End 

Substructure 


Start 

Substructure 


/ 10l(l)l02 (4) 105 


' Z r 1 03(3)l04 / 


y 0 1 ” A. 

/ / (9) 109 110 


111 @ | \® ”6 

- __ 1 ■ 3 j \ 


( 5 }, ^0) Curved Plate Elements 
(i) General Beam Element 

(l2) Rectangular Beam Element 

(l6) Circular Laminated Beam Element 


Rest are Flat Plate Elements 


Figure 8. 1 Example No. 1 ■ Potpourri Input Problem 



POTPOURRI INPUT PROBLEM 




1 

2 


17 

3 

13 2 0 


500. 

6 

2 

11 0 

1 


N 101 

1 

0.0 

Q<0 


N 102 


1.0 

0.0 


H 103 


1.0 

2.0 


N 104 


2.0 

2.0 


H 103 


5.0 

0.0 


N 106 


15.0 

0.0 


N 107 

4 

21.0 

-0.7 


H 108 


21.0 

0.0 


N 109 


30.0 

0.0 


N 110 


36.0 

0.0 


N HI 


47.0 

0.0 


N 112 


60.0 

0.0 


H 113 

3 

60.0 

14.0 


N 114 


66.0 

0.0 


N 115 


92.0 

0.0 


N 116 


92.0 

1.0 


N 117 

2 

1X0.0 

1.0 


D 3 

101 

106 107 115 

114 117 


3 107 


1. ED5 



E 115 

116 

108 100 111 

111 114 114 


C 105 
F 6 

112 

0.0 0.0 

2.0 

0.0 

F 7 


0.0 0.0 

0.0 

2-0 

F 9 


2-0 0.0 

0.0 

0.0 

F 18 


-1-0 

0.0 


THICK 


.02.. 04, 0.05 / 


TABLE 

3 

30. £6 (30.ES , .3 ■ 

11.5E6/ 


TABLE 

2 

10.EB.9.8C6..3 

#3 .805 / 


TABLE 

1 

15.ES.15.E6..3 

, 5.0E6 / 


F 1 

101 

102 1 

1 



.1 

1.E7 1 

.E7 .3 

3.E7 


-1.0 -1.2 -1.3 

-.3 -.25 -.23 

P 2 102 103 1 3 1 

<1.3 .2/3. 1.2) 

-1.3 -1.2 -1.1 

-.25 -.27 -.3 

P 3 103 104 1 1 

•1 1.E7 1.E7 .3 3.E7 

-.5 0.2 0.5 

-.1 0.0 .1 

P 4 102 105 1 1 

• 1 1.E7 l.CT .3 3.E7 

-.» -.7 -.3 

-.23 -.2 -.1 

P 5 105 10$ 2 0 2 

•1 1.C7 1.E7 .3 3.C7 


0.0 0.0 


0 13 


0 13 


0 13 


0 13 


0 13 -« 0 . 
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3 


t2. c» 


.03 

30. E 6 

-.2 

-.2 -21 

-.05 

-.03 -.06 

P 6 

106 106 1 1 2 

(3 .3/3 

.1) 

0.0 

.03 .07 

-.2 

-.1 -.3 

.1 

.02 -03 

-.1 

-.25 -.3 

0.1 

.12 .13 

0.0 

-0.03 -.02 

P 7 

107 108 1 1 1 

.1 1. E 7 1. E 7 

-.1 

-.15 -.15 

-.01 

-.02 -.015 

-.1 

-.11 -.12 

0.0 

.01 -02 

-.2 

-.3 -.5 

.01 

.02 .03 

p 9 

108 109 111 

.1 1 . E 7 1 « E 7 

-.27 

-.29 

-.1 

-.1 

-.3 

-.23 

-.05 --05 

-.36 


0.0 

0.0 

P 10 

109 110 2 0 1 

.1 1. E 7 1. E 7 

-.33 

-uo 

-.2 

-.3 

-.3 

-.2 

-.3 

-.2 

-.3 

-.2 

-.3 

-.2 

P 11 

110 111 1 1 
.1 1 - E 7 1. E 7 


-.6 

-.1 

-.3 

-.5 

-.3 

-.1 

-.12 

-.3 

-.5 

-.1 

-.12 

P 13 

111 112 1 1 
.1 1. C 7 1. E 7 

-.3 

-.6 


-.23 


0.0 

14 


.23 


66 

-.65 

33 

-.1 

66 

-.65 


112 


.1 

113 


3. E 7 


3. E 7 


3. E 7 


.3 


3. E 7 


.3 


3 .ET 


0 2 3 


0 2 3 


0 2 


0 2 13 • 


0 2 


0 2 


0 2 
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1 


l.CT 


l.cr 


3.E7 


-.a -.9 

-.1 -.2 

-.85 -.60 

-.65 -.60 

Q.O 0.0 

0.0 .1 

P 15 112 114 1 2 

Q, 3/3.1) 

-.25 -.395 

-.01 -.03 

P l7 114 115 1 1 

.1 1.E7 1.E7 

-.95 -1.35 

-.3 -.3 

P 18 116 117 1 1 2 

.1 1.E7 1.C7 

.3 1.E7 1.E7 

-1.0 -1.2 -1-3 

-.1 -.2 -.3 

B 8 108 1 1 

1 .67 3.E6 .1 

- 10. 0 
- 6.0 
-7.0 

B 12 111 2 1 1 

.2 1.E7 ‘ 3.E7 

-6.5 


-1.0 



-.3 

16 114 

3 3 

1 

0.2 

1.E7 

3 .06 

0.23 

15. E7 

3.5 E6 

0.3 

30. E6 

12. E6 


.3 3.E7 


.3 3.E7 

.3 3.E7 

-1.5 
-.5 


0 3 


0 3 


0 3 5 


0.5 

.15 2.0 0.0 
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8.2 Example 2: Thermal Stress Analysis of a Stiffened Plate 

The geometry and the temperature distribution of a stiffened 
plate is shown in Figure 8.2. In the cross sectional view, the 
node numbers are indicated as uncircled numbers and the element 
numbers are circled. The local axes of each plate element is 
shown with the subscripts indicating the element. The ends and 
sides of the panel are all simply supported. The properties of 
the plates and beams are shown below. 

Plate Elements 

Thickness = 0.25 in. 

E = 1 x 10 7 lbs/in 2 

v = 0o3 

a = 1 x 10"^ in/in/V 

Beam Elements 

2 

Area = 1.0 in 

E = 1 x 10 7 lbs/in 2 

G - 3.85 x 10 6 lbs/in 2 

a = 1 .0 x IQ' 5 in/in/^ 

The temperature distribution of the panel along its width is 
parabolic with the peak temperature of 370°F at the middle of 
the panel and side temperatures of 10°F o As discussed previously, 
there is no variation of temperature in the Xg or x directions. 
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The coefficients of the polynomial, Tp^ = a y + + c^, that 

describes the temperature distribution for each plate element 
are listed below. 

Plate Element 



a 

b 

c , 


_X_ 

_X_ 


1 

210.0 

13.333 

-.27778 

3 

360. 

3.333 

-.27778 

4 

360. 

-3.333 

-.27778 

6 

210. 

-13.333 

-.27778 


The beam elements, idealized with axial stiffnesses only, are 
input as General Beam elements* The required thermal axial 
load (EaAT) for the beam subjected to AT of 330°F is calculated 
from Equation (4.53) of Ref. 1 as 33000 lbs. 

The harmonics of the Fourier series expansion for the constant 
temperature in the axial direction is truncated at 31. 

If desired, only half of the panel along the width could have been 
analyzed due to symmetry. 
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EXAMPLE NO-2 STIFFENED PLATE 


a o 

0 

0 

0 0 

0 


0 0 



2 

0 29 

31 10 

5 2 

4 

0 



48.0 

4 2 

1 

1 




N 1 1 



0 . 

0 . 


N 2 0 



24. 

0 . 


N 3 0 



36. 

0 . 


N 4 0 



48. 

0. 


N 5 1 



72. 

0 . 


E 2 2 

4 

4 




D 1 1 

5 





PI 1 

2 

1 

0 1 

-0 

-0 -0 -0 

.25 

1 . 

E? 

1 . ET 

.3 

3.8462E6 

1 .0 

210. 

13 

.333333 - 

.2777777 

1 .0 

1. E-5 

1 . 

E-5 




P 3 2 

3 

1 

0 1 

-0 

-0 -0 -0 

.25 

1 . 

ET 

1. ET 

.3 

3.8462E6 

1.0 

360. 

3. 

3333333 - 

.2777777 

1.0 

1. E-5 

1 . 

E-5 




P 4 3 

4 

1 

0 1 

-0 

-0 0-0 

.25 

1 . 

ET 

1. ET 

.3 

3.8462E6 

1.0 

360 . 

-3 

.333333 - 

.2777777 

1.0 

1. E-5 

1 . 

E-5 




Pi 4 

5 

1 

0 1 

-0 

-0 0-0 

.25 

1 . 

ET 

1 . ET 

.3 

3.8 4 62 E6 

1.0 

210. 

-13.33333 - 

.2777777 

1 .0 

1. E-5 

1. 

E-5 




B 2 2 

1 

-0 

0 

1 1 

1.0 

1. ET 

3.85 

E6 

0.0 

0.0 





33000. 



B 9 4 

1 

-O 

0 

1 1 

1.0 

1. ET 

3.85 

E6 

0.0 

0.0 





33000. 




1 1 


1 1 


1 1 


1 1 


31 
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PROGRAM S 0 3 4 4 A / BUCUASP3 


NOV 01 72 


THERMAL STRESS and BUCKLING of LA Ml M«.T£B stiffened panels 
LOADING — IN-FLANE BIAXIAL , THERMAL 

BOUNDARY CONDITIONS — (1) SIMPLY SUPPORTED ALONG EDGES PARALLEL TO THE TRANSVERSE AXIS 

(2) ELASTICALLY RESTRAINED ALONG ANY EDGE PARALLEL TO THE LONGITUDINAL AXIS 
TYPES OF ELEMENTS — (1) PLAT PLATE 

12) CURVED PLATE 

(3) BEAM 

EXAMPLE NO. 2 stiffened plate 


THIS IS AN EDITTED EXAMPLE OF THE PROGRAM OUTPUT. 
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PROGRAM control options 


THERMAL STRESS ANALYSIS 


PROBLEM CHARACTERISTICS 

NUMBER CP NODES = 5 

NUMBER OF FLAT PLATE ELEMENTS = 4 

NUMBER CP CLPvEC PLATE ELEMENTS = 0 

NPJMBER CP BEAM ELEMENTS = 2 

PANEL LENGTH = 46.000 

MAXIMUM PCURIER HARMONIC = 31 

ANALYSIS OPTIONS 

LOAD OPTION (LCACCP) = 4 

WAVE NUMBER SEARCH OPTION (MOPT) = 2 

INITIAL VALUE OF LONGITUDINAL WAVE NUMBER (MMl) = 1 

FI ML VALUE OP LONGITUDINAL WAVE NUMBER (MMA ) = 1 

LIMIT FOR ROOT SEARCHING CRITERIA (NLW) = 0 

UFFER LIMIT FOR ROOT SEARCHING CRITERIA (NUF) = 1 
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FLAT PLATE ELEMENT USER NUMBER 1 USER NOCE I 

NUMEER OF BUCKLING SUBELEMENTS = 2 

NUMBER OF RELATIVE DISPLACEMENT SUBDIVISIONS - 5 

material PROPERTY INFUT OPTION = -0 

MATERIAL PROPERTY I NPL|T FCRMAT OPTION = -0 


USER NOCE J 


THERMAL DATA INPUT OPTION 
LAYER thickness 

1 2.5O0a0OE-Ql 

A- MATRIX 

2747252.747 

624175.624 

0.000 


B~ MATH I X 


D“ MATRIX 


0.000 

0.000 

0.000 


14306.606 

4292.562 

0.000 


EXX EYY 

1.0G0GGGE+07 * 1.000Q00E+07 


824175.824 
274 7252.747 
0.000 


0.000 
0.000 
0 .000 


4292.582 

14308.608 

0.000 


NUXY 

3.0G0000E-01 


MJY X 

3.00GGGOE-01 


0.000 

0.000 

961550.000 


0.000 

0.000 

0.000 


0.000 

0.000 

5008.073 


THERMAL PARAMETERS 


TX 


TTA 


TYB 


TYC 


TZA 


l.OOOCOE+OO 2.100G0E+02 1 .33333E+01 -2.77778E-01 1 .OQOGOE+OG -0 . 

later alfhx alfhy rz 

1 1.00000E-05 0. -0. 


TZB 


8J2 


G 

3.846200E+06 


TZC 


I 



* 

THERMAL * 

STRESSES AND DISPLACEMENTS * 

* 

*44 4 $ 4 $ £ 4 4 $ $ $ $ $ 4 * 44 44 * 44 44 4 4 $ $ 4 4 44 44-4 44444 * 44 44444 4444 4 * 4 44 444 4 


FOURIER TERM NUMBER = 29 


START SUBSTRUCTURE NUMBER - 


PLATE ELEMENT NUMBER - 1 


X 

U 

V 


w 

Nil 

F£2 

M2 

Ml 


0.0000 

7.3S19E-15 

0. 

0. 


0. 

0. 

-8.3078E+03 

0. 

0. 

2. 6567 

7.2363E-15 

-1 .0540E-G2 

0. 


-2.3985E+02 

9 . 6492E-1 1 

-6.8467E+03 

0. 

0. 

3.3333 

6.8946EM5 

-1 .9221E-02 

0. 


-2.58$0e*Q2 

2.3493E-10 

-3.71Q5E+03 

0. 

0. 

8.0000 

6.321 7E-1 5 

-2.6393E-02 

0. 


-2.6053E*02 

3 . 7949 E- 10 

-4.7225E*03 

0. 

0. 

10.6667 

5.59SOE-15 

-3.2233E-02 

0. 


-2.5383E+02 

4. 3119 E - 10 

-3.8216E+03 

0. 

0. 

13.3333 

4.720 7E-15 

-3.6865E-02 

0. 


-2.4 53SE+Q2 

5.4629E-10 

-2.978GE+03 

0. 

0. 

16.0000 

3.6481E-15 

-4.0381E-02 

0. 


-2.4389E+02 

6.53S7E-10 

-2.179SE+03 

0. 

0. 

18.6667 

2.4923E-15 

-4.2849E-02 

0. 


-2.4724E+02 

7.1245E-10 

-1 .4226E+03 

0. 

0. 

21 .3333 

1 . 2 708 E' 1 5 

-4.4314E-02 

0. 


-2.5272E+02 

7.3342E-10 

-7.0075E+G2 

0. 

0. 

24.0000 

9.0656E-29 

-4 .4799 E- 02 

0. 


-2.5S30E+02 

7. 4090 E- 10 

-4.4614E-H 

0. 

0. 

QURIER TERM NUMBER = 

29 START 

SUBSTRUCTURE 

NUMBER = 1 

FLATE ELEMENT NUMBER 

= 1 

Y = 

X 

U 

V 


u 

Ml 

N22 

M2 

Mil 


0.0000 

-2.4493E-02 

0. 

0. 


0. 

0. 

-7.3659E+C3 

0. 

0. 

2.6667 

-2.1026E-02 

-9.5791E-03 

0. 


-3.e096£>02 

-2.4965E+03 

-6.40C2E 4 03 

0. 

0. 

5.3333 

-1 .76S6E-02 

-1 .771 6E-02 

0. 


-6.7223E+G2 

-2.43O8E*03 

-5.4090E+03 

0. 

0. 

8.0000 

-1 .4631 E-02 

-2.4512E-02 

0. 


-8.9304E + 02 

-2.2271E+03 

-4.4929E+G3 

0. 

0. 

10.6667 

-1 .1854E-02 

-3.0072E-02 

0. 


-1 .0613E+03 

-1 .9 731 E-H33 

-3.6423E+D3 

0. 

D. 

13.3333 

-9.2576E-03 

-3.4494E-G2 

0. 


-1 .1635E+03 

-1.75&3E+G3 

-2.8445E*03 

0. 

0. 

16.0000 

-6.801 7E-03 

-3.78 54E-02 

0. 


-1 .2818E+03 

-1 .6211 E+03 

-2.0898E+G3 

0. 

0. 

18.6667 

-4.4603E-03 

-4.0214E-02 

0. 


-1.3457E^03 

-1 .5925E+03 

-1 .3705E+G3 

0. 

0. 

21.3333 

-2.2059E-03 

-4.1614E-02 

0. 


-1.383GE+03 

-1 .6177E+03 

-6.7766E+02 

a. 

0. 

24.0000 

-1.3035E-16 

-4.2078E-02 

0. 


-1.3952E*G3 

*1 .6353E+03 

-4.0075E-11 

0. 

0. 


- 12.0000 

K22 
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V 


thermal stresses and displacements <CC#JT) 


FCUUER term number = 

29 START 

SUB STRUCT LIRE 

NUMBER = 1 

FLATE ELEMENT NUMBER 

= 1 

Y = 

0.0000 

X 

U 

V 

W 

Nil 

^2 

N3 2 

Mil 


K22 

0.0000 

-4.5992E-02 

0. 

0 . 

0 . 

0 . 

-6.271 1 E + 03 

0. 

0. 


3.6667 

-3.9811 E~02 

-0.1756E-O3 

0. 

-5.1343E+02 

-4 .5300E+03 

-5.4919E+03 

0. 

0 . 


3.3333 

-3.3633E-02 

-} .5170E-02 

0. 

-9.4194E+02 

-4.4993E+03 

-4 . 6775E+03 

0. 

0. 


8 .0000 

-2.7981E-02 

-2.1042E-02 

0. 

-1 .2938E+03 

-4.1747E+03 

-3.9lG6E + 03 

0 . 

0. 


10.6667 

-2.2736E-02 

-2.5864E-02 

0. 

-1 . 5750E+03 

-3.7371E+03 

-3.1878E+03 

0. 

0. 


13.3333 

-1 .7797E-G2 

-2.9710E-02 

0. 

-1.7924E+03 

-3.3453E*03 

-2.3015E+03 

0 . 

0. 


16.0000 

-1 .3101E-02 

-3.2 639E-02 

Q. 

-1.9530E+03 

-3.11 53E 4 03 

-1 .6452E+03 

0. 

0. 


18.6667 

-8.6048E-03 

-3.4699E-02 

0. 

-2.0631E+03 

-3.062GE+03 

-1 .21 39E 4 03 

0. 

0. 


21 .3333 

-4 .2601 £-03 

-3.5923E-02 

0. 

-2.1274E+03 

-3 .1G42E+03 

-6.0163E+02 

0. 

0. 


24.0000 

•2 * 5398 £-1 6 

- 3. 63 29 E- 02 

0. 

-2.1486E+03 

-3.1343E+03 

-3.6226E-11 

0. 

0. 



FCXRIER TERM NLP-DER - 

29 START 

SUBSTRUCTURE 

NUMBER = 1 

PLATE ELEMENT NUMBER 

= 1 

r = 

6.0000 

X 

U 

V 

W 

Nil 

N22 

N12 

Mil 


M22 

o.oooo 

-6.3593E-02 

0 . 

0 . 

0. 

0 . 

-4.9202E+03 

0 . 

0 . 


2.6667 

-5.S297E-02 

-6.4718E-03 

0 . 

-6.2683E+02 

-6.1 566E+03 

-4.31 S7E+03 

0 . 

0 . 


5.3333 

-4.€935E~02 

-1 .19/GE-02 

0 . 

-1 .1648E+03 

-6.1556E+03 

-3.691GE+G3 

0 . 

0 . 


8.0000 

-3.9226E-Q2 

-1 .6562E-02 

0 . 

-1 .5926E+03 

-5.7384E+03 

-3.1162E+03 

0 . 

0 . 


10.6667 

-3 .2011 E-02 

-2.0 338E-02 

0 . 

-1 .9208E+03 

-5.1 713E f 03 

-2.569GE+03 

0 . 

0 . 


13.3333 

-2.5149E-02 

-2.3361E-02 

0 . 

-2.1696E+03 

-4.6655E+03 

-2.035SE+03 

0 . 

0 . 


16.0000 

-1 .856SE-02 

-2.5675E-02 

0 . 

-2.3532E+03 

-4.3725E+03 

-1 .5127E+03 

0 . 

0 . 


18.6667 

-1 .2221 E“ 02 

-2 . 7308 E* 02 

0 . 

-2 .4800 E* 03 

-4 .31 11 E + G3 

-1 .00036+03 

0 . 

0 . 


21.3333 

-6.05S6CHj3 

-2 .8260E-02 

0 . 

-2.5546E+03 

-4.3727E+03 

-4.9 724E+02 

0 . 

0 . 


24.0000 

-3.6543E-16 

-2.8603E-02 

0 . 

-2.5792E+03 

-4 .4146E+03 

-3.0464E-11 

0 . 

Q. 



rOLRIER TERM m>£ER = 

29 START 

SUBSTRUCTURE 

NUMBER = 1 

PLATE ELEMENT NUMBER 

= 1 

Y = 12 .0000 

X 

U 

V 

W 

Nil 

NC2 

N12 

Mil 

M22 

0.0000 

-7.5061E-02 

0. 

0. 

0. 

0. 

-1 .9861 £+03 

0. 

0. 

2.6667 

-6.6569E-02 

-4.3467E-03 

0. 

-1.4582E+03 

-7.0755E t 03 

-3.1428E+03 

0. 

0. 

3.3333 

-5.7272E-C2 

-8.1034E-03 

0. 

-1.7170E+03 

-7.1889E+G3 

-2. 82 79 C + 03 

0. 

0. 

6.0000 

-4.828 7E-02 

-1 .1293E-02 

0. 

-1.9501E+03 

-6.6085E+G3 

-2.38 79E+03 

0. 

0. 

10.6667 

-3.9S31E-02 

-1 .394GE-02 

0. 

-2.1265E+03 

-6.2271E+03 

-1 .9507E+03 

0. 

0. 

13.3333 

-3.1264E-02 

-1 . 60 70 r- 02 

0. 

-2.2736E«03 

-5.6895E+03 

-1.5489E+03 

0. 

0. 

16.0000 

-2.3156E-02 

-1 .77O2E-02 

0. 

-2.4229E*03 

-5.3751 E*03 

-1 .1739E+03 

0. 

0. 

18.6667 

-1 .527SE-C2 

-1 .8855E-02 

0. 

-2 . 5746E 4 03 

-5.311 €E 4 03 

-8.0023E+G2 

0. 

0. 

21-3333 

-7.S867E-03 

-1.9542E-02 

0. 

-2.6949E*03 

-5.3821E*03 

-4.O068E+O2 

0. 

0. 

2 4.0000 

~4. 6441 E- 16 

-1 .9770E-02 

0. 

-2.7413E*03 

-5.4292E+03 

-7.3702E-12 

0. 

0. 
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THERMAL STKCSSCS A HI C I S PLACEMENTS (COW) 


FOURIER TERM fWJMDER = 

29 START 

SUBSTRUCTURE 

II 

«S 

9 

2 

beam element number = t 

X 

MZ 

MT 

TCfctfUE 

AXIAL load 


0.0000 

0. 

o'. 

0. 

0. 


2.666? 

0. 

0. 

0. 

2 . 6578E 4 03 


5.3333 

0. 

0. 

0. 

1.7589E+03 


e.oooo 

0. 

0. 

0. 

3.6974E+02 


10.666? 

0. 

0. 

0. 

-1 .0336E+03 


13.3333 

0. 

0. 

0. 

-2.2670E+03 


16.0GGG 

0. 

0. 

0. 

-3.2414E+03 


18.666? 

0. 

0. 

0. 

-3.9243E+03 


21 .3333 

0. 

0. 

0. 

-4.3209E+03 


24.0000 

0. 

0. 

0. 

“4.4500E + 03 



FOURIER TERM NLRGER = 

29 START 

SUBSTRUCTURE 

NUMBER = 1 

PLATE ELEMENT NUMBER 

= 3 

T = -6.0000 

X 

U 

V 

U 

Nil 

N£2 

N12 

mi 

K22 

O.OGGO 

-7.50 61E-02 

0. 

0 . 

0 . 

0 . 

-4.66G9E+03 

0 . 

0 . 

2.6657 

-6.6569E-02 

-4 .3467E-03 

0 . 

-1 .4582E+03 

-7 ,0755E 4 03 

-2.7229E+03 

0 . 

0 . 

3.3333 

-5.72 72E-02 

-8. 1034E-03 

0 . 

-1 .7170E*03 

-7.1889E+03 

-2.2130E+03 

0 . 

0 . 

8.0000 

-4 .8267E-02 

-1 .129 3E-02 

0 . 

-1 .9501E+O3 

-6.8O8 5E+03 

-1 .8512E+03 

0 . 

0 . 

10.6667 

-3.9631 E-02 

-1 .3940E-02 

0 . 

-2.1265E+03 

-6.2271 E + G3 

-1 . 5355E*03 

0 . 

0 . 

13.3333 

-3.1264E-02 

-1 .5070E-02 

0 . 

-2. 2 73 6E 4 03 

-5.6S95E + 03 

-l .1996E + 03 

0. 

0 . 

16-OOGO 

-2.315SE-C2 

-1 .7702E-O2 

0 . 

-2.4229E+03 

-5-3751E+03 

-8.6004E+02 

0 . 

0 . 

18.6667 

-1.5279E-02 

-1 .88 55E-02 

0 . 

-2. 5?4SE + 03 

-5.3116E + 03 

-5.4065E+02 

0 . 

0 . 

21.3333 

-7.58 67^-03 

-1 .9 542E-02 

0 . 

-2.6949E*C3 

- 5.3821 E+03 

-2.56S9E+02 

0 . 

0 . 

24.0000 

-4 .6-44 IE- 1 6 

-1 .9770 E'02 

0 . 

-2.7413E+03 

-5.4292E+03 

-3.3551 E~ 11 

0 . 

0 . 


FOLRIER TERM NUT'S er = 

29 START 

SUBSTRUCTURE 

fJLTCER = 1 

PLATE ELEMENT NUMBER 

= 3 

T 2 -3.0000 

X 

U 

V 

W 

Nil 

t£2 

N12 

mi 

►C2 

0.0000 

-6.0982E-02 

0. 

0. 

0. 

0. 

-2.5117E+03 

0. 

0. 

2.6667 

-7.0S95E-02 

-3.112EE-03 

0. 

-8.3005E+02 

-7.7232E+03 

-2.2 4 62E+03 

0. 

0. 

5.3333 

-6.0328E-02 

-5.9842E-03 

0. 

-1.4397E+03 

-7 . 7285E*03 

-1 .87220+03 

0. 

0. 

6.0000 

-5.0665E-02 

-8.4718E-03 

0. 

-1 .8 658E+D3 

-7.2590E+03 

-1 .519GE+03 

0. 

0. 

10.6667 

-4.1493E-02 

-1 .0541E-02 

0. 

-2.195ZE*03 

-6.6012E+03 

-1 .2101 E + 03 

0. 

0. 

13.3333 

-3.2 684E-D2 

-1 .2201E-02 

0. 

-2.4 5 63E+03 

-6.0047E+03 

-9.3433E+02 

0. 

0. 

16.0000 

-2.4176E-02 

-1 .3471 £-02 

0. 

-2.6570E+03 

-5.6EG0E+03 

-6.81 55E+02 

0. 

0. 

18.6667 

-1 .5933E-02 

-1 .43660*02 

0. 

-2.7997E+03 

- 5 . 5944 E*03 

-4.4512E+02 

c. 

0. 

21 .3333 

-7.9050E-03 

-1 .48 98E-02 

0. 

-2.S850E+03 

-S.6763E*03 

-2.1971 E + 02 

0. 

0. 

24.0000 

-4.8335E-16 

-1 .5075E-02 

0. 

-2.9134E*03 

-5.7300E*03 

-1 .15DGE-11 

0. 

0. 
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8.3 Example 3: Thermal Buckling of a Cylinder Heated Along an Axial Scrip 

Figure 8.3 shows the cylinder properties and the idealized structure 
of this sample problem which also appears in Ref. 3. The problem 
here is to compute the critical temperature ratio, the value of Tq 
when buckling occurs. 

The circled numbers indicate element numbers and the node numbers 
are uncircled. Only half of the cylinder is considered with 
appropriate symmetry conditions imposed. As the figure indicates, 
three groups of elements are used, all elements equally spaced 
within its group. Every element is divided into 5 subelements. 

The initial temperature distribution in the hoop direction is 
approximated by seven different parabolic equations for elements 
(T) to (?) and no temperature change for the rest of the elements. 


8.16 






EXAMPLE MO. 3 STRIP HEATED CYLINDER 



0 

0 

0 

0 0 0 0 


0 

1 

1 

2 5 49 


21 

0 

0 

20 


6 

2 

139 

139 

N 

1 

4 


0.0 5.181 

N 

2 

0 


.18081 <29 5 . 1 778439 

N 

3 

0 


.361406295.1663793 

N 

4 

0 


.541561975.1526179 

N 

5 

0 


.7210 5 56 4 5.130 5789 

N 

6 

0 


1.0771905 5.0677827 

N 

7 

0 


1.4260771 4.9602968 

N 

8 

0 


1 .7720064 4.8685475 

N 

9 

0 


2.769057 4.37893 

N 

10 

0 


3.638845 3.68803 

N 

11 

0 


4.34136 2.827606 

N 

12 

0 


4.844319 1.837208 

N 

13 

0 


5.124604 0.762361 

N 

14 

0 


5.169331 -0.347529 

N 

15 

0 


4.976444 -1.441445 

N 

16 

0 


4.55481 -2.469103 

N 

17 

0 


3.923809 -3.383265 

N 

18 

0 


3.112446 -4.141913 

N 

19 

0 


2.156016 -4.710173 

N 

20 

0 


1.104390 -5.061925 

N 

21 

4 


0.0 -5.181 

C 

2 

2 

10 

11 20 

S 

1 



1. £20 1. £20 

s 

21 



1. £20 1. £20 

p 

1 

1 

2 

2 0 1 -O 


3. E 7 .3 

■ .3258147 -1 . 755601 


,006 3. ET 

1.0 .904664 5 

6. E-6 6. E-6 

2 2 3 2 0 1 H 

.006 3. ET 3. E7 .3 

1.0 .8719024 -.6S9G4 55 -1.172957 

6. E-6 6. E-6 

3 3 4 2 0 1 H 

.008 3. E7 3. E7 .3 

1.0 .6833338 -1.137967 -.3066789 

6. E-6 6. E-6 

4 4 5 2 0 1 H 

.008 3. E7 3. E7 

1.0 .4741115 -1.109465 

6. E-6 6. E-6 

5 5 6 2 0 1 

.008 3. E7 3. E7 

1.0 .2180377 -.7346457 

6. E-6 6. E-6 

6 6 7 2 0 1 

.008 3. E7 3. E7 

1.0 .05052683 -.2524012 


51 

48.0 


.3 

.426744 

-0 

.3 

.6064855 

-0 

.3 

.4702425 


-O 0 
11 .53E6 
1.0 

-O 0 
11.53E6 
1 .0 

-O 0 
11.53ES 
1.0 

-0 0 
11.53E6 
1.0 

-0 0 
11 .53E6 
1.0 

-0 0 
11.53E6 
1 .0 


-0 


-0 


-0 


51 


-5.181 


-5.181 


-5.181 


-5.181 


-5.181 


-5.181 


8J8 



6. E-6 6. E-6 


p 

T T 

8 2 

0 

1 

-0 

i 

o 

o 

-0 

1 

1 

-5.181 


.008 

3. ET 

3. 

E7 

.3 

11 .53E6 






1 .0 

007191883- 

.04976467 

.1243829 

1 .0 






6. E-6 

6. E-6 









p 

8 6 

9 2 

0 

1 

-0 

-a 0 

-0 

1 

1 

-5.181 


.008 

3. E7 

3. 

E7 

.3 

11 .53E6 






1 .0 

1 .0 










6. E-6 

6. E-6 









p 

9 9 

10 2 

0 

1 

-0 

-0 0 

-0 

1 

1 

-5.181 


.008 

3. E? 

3. 

E7 

.3 

11 .53E6 






1.0 

1.0 










6. E-6 

6. E-6 









p 

10 10 

11 2 

0 

1 

-0 

* 0 0 

-Q 

1 

3 

-5.181 


.008 

3. E7 

3. 

E7 

.3 

1 * . 53Et> 






1.0 

1.0 










6. E-6 

6. E-6 









p 

11 11 

12 2 

0 

1 

-0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. E7 

3. 

ET 

.3 

11 .53E6 






1.0 

1.0 










6. E-6 

6. E-6 









p 

12 12 

13 2 

0 

1 

-0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. ET 

3. 

ET 

.3 

11.53E6 






1.0 

1.0 










6. E-6 

6. E-6 









p 

13 13 

14 2 

0 

1 

“0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. ET 

3. 

E7 

.3 

11 .53E6 






1.0 

1.0 










6. E-6 

6. E-6 









p 

14 14 

15 2 

0 

1 

-0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. ET 

3. 

ET 

.3 

11.53E6 






1.0 

1 .0 










6. E-6 

6. E-6 









p 

15 15 

16 2 

0 

1 

-0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. ET 

3. 

ET 

.3 

11.53E6 






1.0 

1.0 










6. E-6 

6. E-6 









p 

16 16 

17 2 

0 

1 

-0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. E7 

3. 

E7 

.3 

11.53E6 






1 .0 

1.0 










6. E-6 

6. E-6 









p 

IT 17 

18 2 

0 

1 

-0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. E7 

3. 

E7 

.3 

11.53E6 






1.0 

1.0 










6. E-6 

6. E-6 









p 

10 18 

19 2 

0 

1 

-0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. E7 

3. 

E7 

.3 

11 .53E6 






1.0 

1.0 










6. E-6 

6. E-6 









p 

19 19 

20 2 

0 

1 

-0 

-0 0 

-0 

1 

3 

-5.181 


.008 

3. E7 

3. 

E7 

.3 

11 .53E6 






1.0 

1.0 
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0 1 
3. zr 


l 


3 


-5.181 


6. E-6 6. E-6 

p 20 20 21 2 

.008 3. E7 

1.0 l.o 

6. E-6 6. E-6 

L 


-0 -o 0 -o 

.3 11.53E6 
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F R O 0 R A M 3 0 3 4 4 A / BUCLASP3 NOV 01 72 

THERMAL STRESS AND SUCKLING OF LA Ml MATED STIFFENED PANELS 

LOADING — 1N-FLANE BIAXIAL , THERMAL 

BOJNCARY CONDITIONS — (1) SIMPLY SUPFORTED ALONG EDGES PARALLEL TQ THE TRANSVERSE AXIS 

(2) ELASTICALLY RESTRAINED ALONG ANY EDGE PARALLEL TO THE LCNGI TUCI NAL AXIS 
TYFES OF ELEMENTS — (1) FEAT PLATE 

(2) CURVED FEATE 

(3) BEAM 

EXAMPLE NO. 3 STRIP HEATED CYLINDER 


THIS IS AN EDITTED EXAMPLE OF THE PROGRAM OUTPUT. 


8.21 



NODAL 


DATA 



COORDINATES 

BOUNDARY CONDI TI CMS 


SPRING STIFFNESSES 

NODE 

Y 

2 

CODE 

MW) 

K (THETA) MV) K 

1 

0 .QGOOQ 

5.18100 

4 

-0. 

1 .OOOGQOOOE+20 1 .OOOOODGOE+20 -0 

2 

.18081 

5.17734 

0 



3 

.36141 

5.16838 

0 



4 

.54156 

5.15262 

0 



3 

.72106 

• 5.130 53 

0 



6 

1 .07719 

5.06778 

0 



7 

1 .42808 

4.98030 

0 



A 

1.77201 

4.86855 

0 



9 

2.76307 

4.37893 

0 



10 

3.6388 5 

3.66303 

0 



11 

4.34136 

2.82761 

0 



12 

4.84432 

1 .8 3 721 

0 



13 

5.12460 

.76236 

0 



14 

5.16933 

-.34753 

0 



15 

4.9 7644 

-1 .44145 

0 



16 

4.55481 

-2.46910 

0 



17 

3.92381 

-3.38326 

0 



18 

3.11245 

-4.14191 

□ 



19 

2.15802 

-4.71017 

0 



20 

1 .10439 

-5.06193 

0 



21 

0 .00000 

-5.18100 

4 

- 0 . 

1 . 00000000 E* 20 1 . 00000000 r+ 20 -0 . 

S U B S T 

R U C T U R E 

D E F I 

NITION DATA 



MJMEER OF 

SUBSTRUCTURES 

= 2 





SUBSTRUCTURE FIRST NODE LAST NjCE 


START 

2 

10 

END 

11 

20 


r 
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ELEMENT TRANSFORMATION DATA 


ELEM. NO. 

WIDTH 

transformation 
SIN(I) COSU) 

SIN(J) 

CCMJ> 

EL. TYPE 

: 

.18085 

.00000 

1 .00000 

-.03490 

.99939 

CRVD 

PLATE 

2 

.18085 

-.03490 

.99939 

- ,06976 

.99756 

CRVD 

PLATE 

3 

.18085 

- .06976 

.99756 

-.10453 

.99452 

CRVD 

FLATE 

4 

.18085 

-.10453 

.99452 

-.13917 

.99027 

CRVD 

FLATE 

5 

.36170 

-.13917 

.99027 

-.20791 

.97815 

CRVD 

FLATE 

6 

.36170 

-.20791 

.97815 

-.27564 

.96126 

CRVD 

FLATE 

7 ■ 

i 

.36170 

-.27564 

.96126 

-.34202 

.93969 

CRVD 

FLATE 

8 

1.11293 

-.34202 

.93969 

-.53447 

.84519 

CRVD 

FLATE 

9 

1.11293 

-.53447 

.84519 

-.70234 

.71184 

CRVD 

PLATE 

10 

1.11293 

-.70234 

.71184 

-.83794 

.54576 

CRVD 

FLATE 

11 

1.11293 

-.8 3 794 

.54576 

-.9 3 502 

.35460 

CRVD 

PLATE 

12 

1.11293 

-.93502 

.35400 

-.98911 

.14715 

CRVD 

FLATE 

13 

1 .11293 

-.98911 

.14715 

-.99775 

-.06708 

CRVD 

FLATE 

14 

1.11293 

-.99775 

-.06708 

-.96052 

-.27822 

CRVD 

FLATE 

13 

1 .11293 

-.96052 

-.27822 

-.87914 

-.47657 

CRVD 

FLATE 

16 

1 .11293 

-.87914 

-.47657 

-.75735 

-.65301 

CRVD 

FLATE 

17 

1 .11293 

-.75735 

-.65301 

- .600 74 

-.79944 

CRVD 

FLATE 

18 

1 .11293 

-.600 74 

-.79944 

-.41652 

-.90912 

CRVD 

FLATE 

19 

1.11293 

-.41652 

-.90912 

-.21316 

-.97702 

CRVD 

FLATE 

20 

1.11293 

-.21316 

-.9 7 702 

.00000 

- 1.00000 

CRVD 

FLATE 

E L E M E N 

T SUB 

S T R U C 

T U R E 

A S S I 

G N M E N T 

D A 1 

' A 


NUMBER BLOCK 

1 START 

2 START 

3 START 

4 START 

5 START 

6 START 
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,002) 7 


. 0000 } 


3 

sirciocs rc* tLC^C‘«i 1 1 i«"iTM,'C t lrc C’i~ 

t -.00290 - .0000} 2 -.00261 -.GCGG3 3 -.00274 -.GGGC3 4 - .00269 - .OOCC 4 

5 -.00265 -.0C0Q4 


nimcss mtkix bancwIct* s c 


»vt>AHC 

STRESSES r & 

CLC>^NT 1 

SUBSTRUCTURE START 


1 

- . 8908 4 

- .00394 

2 -.56424 

-.00392 

3 

-.62722 

- .00370 



AVOiACOt 

STRESSES rcS 

ELEVEN 2 

SUBSTRUCTURE START 


1 

-.79578 

-.0035 a 

2 -.759 76 

-.00345 

s 

-.62 754 

-.00293 



AVE*At€C 

stresses ro< 

ae<vr 3 

4 

SX STRUCTURE START 


a 

-.57375 

- .002 74 

2 -.32264 

-.00253 

s 

-.55622 

-.00188 



AVOiACCE 

stresses res 

D-EXfiT 4 

SUBSTRUCTURE START 


1 

-.5C?t 5 

-.00163 

2 -.24774 

- .00145 

s 

-.09 529 

-.00060 



AVCRA0CC 

STS CSSC3 FC S 

D-CXENTT 5 

SUBSTRUCTURE START 


a 

- .02^4 

- J0GO31 

2 .05329 

-.00016 

9 

•22470 

.000 02 



AVCPAOCC 

STS CSS£3 FCS 

ELCC77T 6 

SUBSTRUCTURE START 


1 

.23630 

.000 78 

2 .27907 

.00091 

9 

.50 55 7 

.00106 



A,V£RAfcCX 

STS CSSC5 FO< 

ELEXLNT 7 

SUBSTRUCTURE START 


1 

.50156 

J00109 

2 . 29609 

.00109 

9 

2690 

.00 1 04 



AVOtAVCT 

STRESSES F0J 

CLOCKT 8 

SUBSTKuCTUvC START 


a 

•24 791 

.OOC-98 

2 .21142 

.00090 

9 

.11659 

.00009 



AVCRAtCE 

sts esses Fes 

ELE*CWT 9 

SUBSTRUCTURE START 


1 

.06969 

.00062 

2 .06533 

.000 55 

9 

•0061 9 

.00032 



ArCRAtCB 

sts esses re* 

ELCXCNT 10 

suesrsuc tiX E start 


1 

- .OG599 

.00023 

2 -.01782 

.00016 

9 

-.0*022 

-.OCOOl 



A rC* A tCE 

s r*. esses fcr 

CLC*EKT 1 

SUBSTRUCTURE CNC 


a 

-.00399 

.00025 

2 -.01762 

.00016 

9 

-.04022 

-.00GG1 



ArCKACOC 

s tr esses FCA 

CLC^N-r 2 

SUBSTRUCTURE ENt 


s 

- .0<379 

-.00006 

2 -.04367 

- .00011 

9 

- .04301 

- .00020 



<irtA< tOC 

stresses ref 

CLC*GMT 5 

substructure CMC 


1 

- .04000 

-.0C022 

2 -.03625 

-.00023 

9 

-.02239 

- .GCG22 



A VCA A CC t 

stresses fca 

cte*<CHT 4 

SUBSTRUCTURE CPO 



3 

-.87143 

-.0038 7 

4 

-.8 32 4 4 

-.00379 

3 

71973 

-.00329 

4 

-.67366 

- .00312 

3 

-.46371 

-.00232 

4 

- . 41 3S 1 

- X0210 

3 

-.19314 

-.00122 

4 

-.14433 

- .00101 

3 

.12413 

.00013 

4 

.18059 

.000*1 

3 

.29480 

,ocr*99 

4 

.33325 

.0910; 

3 

.26893 

.00108 

4 

.26010 

.00106 

3 

.17729 

.00083 

4 

.14560 

.00076 

3 

•04390 

.00047 

4 

.02460 

.00040 

3 

-.02738 

.00011 

4 

-.03481 

.cocos 

3 

-.02738 

.00011 

4 

-.03481 

.OOOOS 

3 

-.04604 

-.00014 

4 

- .04509 

- .GGGt 8 

3 

-.03193 

-.00023 

4 

-.02727 

- .00023 
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1 


t 

-.01746 

-.00021 

2 -.01262 

-.00019 

3 

-.00 799 

-.0C017 

4 

-.00369 

-.00013 

9 

.00023 

-.00012 









*Vt**C,CC 

t 

$TfcC33CS r<* 
.00369 

ri.EK.NT 3 

-.00010 

SUSSTRUCTeKE ENC 
2 <00662 

-.0000? 

3 

.00904 

-.00003 

4 

.01093 

-.00002 

9 

.01229 

.00000 









avcaaocc 

i 

STRESSES r<* 
.01314 

CLEMENT 6 

.00002 

SUBSTRUCTURE ENC 
2 .01332 

,00004 

3 

XI 349 

.00003 

4 

.01309 

.0000. 6 

5 

.0123? 

.00007 









AVEKACCC 

t 

»T*CS3CS r^S 
.01141 

ELEKNT ? 

.00009 

SU&STRuCTURE ENt 
2 .01023 

.00009 

3 

.00996 

.00008 

4 

.00759 

.00009 

9 

.00620 

.00009 









AVCSACCD 

1 

3TSCS3C3 EC* 
.00492 

ELEMENT 6 

.COCO? 

SlXSTRoCnSC ENC 
2 .00350 

.00006 

3 

.0022? 

4X3006 

4 

.00115 

,00005 

9 

.00016 

.00004 









A^dUCtC 

1 

stresses res 

-.00069 

CLEMENT 9 

.00003 

SUBSTRUCTURE ENC 
2 - .00140 

,ffyy>2 

3 

- .0019 7 

4X3001 

4 

-.00240 

.00001 

9 

- .002 n 

- .000 CO 









AvtRACCC 

t 

strcsscs res 
-.00291 

o. CMC NT 10 
- .0000 1 

SUCSTRUCTURE ENC 
2 -.00301 

-.00001 

3 

-X0303 

- .00002 

4 

-430303 

-.00002 

9 

- .0029? 

-X00G3 









AVEAACCC 

1 

s misses res 
- .00290 

CL EMC NT 11 
- .0000 3 

SUBSTRUCTURE ENC 
2 - .00291 

-.00003 

3 

-.00274 

-4X3003 

4 

- .00269 

-.00004 


s -JJMW -_OCCO* 
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EXAMPLE NO-3 STRIP HEATED CYLINDER 


variable ratio cm 

NUMBER Cf NEG. 

BUCKLING DETERMINANT (A* 

FIRST PLATE CLEMENT 

«jN CIAGCMAL 

A 

B 

16G.0000000 

0 

.41021 

1156 

240.000 0000 

0 

.14562 

1156 

260.0000000 

0 

.30644 

1152 

SOO. 000 0000 

1 

-.63932 

1152 

290.0000000 

1 

-.17334 

1152 

26S.OOGOO 00 

0 

.06494 

1152 

28 7 . 50000C0 

1 

-.87369 

1148 

28 6.2500000 

0 

.06090 

1148 

286.8750000 

1 

-.39691 

1148 

286.5625000 

1 

-.15811 

1148 

286.40 62 500 

1 

-.61608 

1144 

266.328125 0 

0 

.33605 

1144 

286.3671875 

1 

-.14004 

1144 

286.3476562 

0 

.09900 

1144 

286.3574219 

1 

-.32838 

1140 

28 6.352 5391 

0 

.62 778 

1140 

286.3549805 

0 

.14970 

1140 

286.3562012 

1 

-.08934 

1140 

286.3555908 

0 

.48292 

1136 

286.3558960 

1 

-.47324 

1136 

286.3557434 

0 

.07745 

1132 

286.355S197 

1 

-.23420 

1136 

23 6.3 5 5 78 1 6 

1 

-.11468 

1136 

286.3557625 

1 

-.87872 

1132 

236-3557529 

1 

-.40063 

1132 

286.3557482 

1 

-.161 59 

1132 

286.3557458 

1 

-.67321 

1128 

286.3557446 

0 

.28297 

1128 

286.3557452 

1 

-.19512 

1128 


COMMENTS 
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8.4 Example 4: Mechanical Buckling of a Heated Plate 

The pertinent data for this example is shown in Figure 8.4. 

The critical N-j-j mechanical load is computed for this plate that 
is heated 35°F. Each plate element is divided into 5 subelements. 
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N 

N 

N 

C 

s 

s 

p 


p 


EXAMPLE NO-4 MECHANICAL BUCKLING QE A HEATED PLATE 


0 0 

0 

0 

0 0 

0 1 



2 

3 0 

2 

0 


5 2 

4 

4 


1 4 



0.0 

2 0 



2.5 

3 4 



5.0 

1 1 

3 



1 

1 . 

E5G 

1 . E5Q 

3 

1 . 

E50 

1. E50 

1 1 

2 

1 

0 1 

.1 

1 

. e y 

1. E? 

1.0 

35 

.0 


1. E-5 

1. 

E-5 


2 2 

3 

1 

0 1 

.1 

1 

. E7 

1. E7 

1-0 

35 

.0 


1. E-5 

1. 

E-5 



0 

2 29 31 0 

15-0 

0.0 

0.0 

0.0 

1 . E50 

1. E50 

-o -o 0 -0 

.3 3-85 E6 

1.0 

-o-o 0-0 

.3 3-85 E6 

1 .0 


L 


1 


1 


1 


1 


31 
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8.5 Example 5: Buckling of a Cylinder Under Bending 

The relevant information concerning this example is shown in 
Figure 8.5. The cylinder is idealized as a half cylinder because 
only symmetrical buckling patterns about the diameter are considered. 
Each element is divided into 5 subelements. The problem here is to 
calculate the critical load ratio which by definition is the ratio of 
the critical bending moment to the initially applied bending moment. 
The initial bending moment is approximated by applying different 
axial loads to each subelement as shown in Figure 8„5 a These axial 
loads, Nii , must be a constant value for each subelement. 

The initial bending moment is chosen such that the crown = 1787 
lbs. /in.. Tensile loads are input as positive values and compressive 
loads as negative values. 
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Idealized Structure 


Load Distribution 


Figure 8.5 Cylinder Under Bending 
8.3T 



EXAMPLE 

NO. 5 

BUCKLING OP 

A CYLINDER 

UNDER BENDING 





0 0 

0 

0 0 0 

0 







0 1 


0 

5 0 

0 0 






5 0 

0 

4 


94.25 




0 


8 2 

16 

19 







N 

1 4 


0.0 

30.0 






N 

2 0 


25.960762 

15.0 






N 

3 0 


30. 

0.0 






N 

4 0 


25.960762 

-15.0 






N 

5 4 


0.0 

-30. 






D 

1 1 

5 








S 

1 


1 . E2G 

1 . E20 






' s 

5 


1 . E20 

1 . E20 






p 

1 1 

2 

2 0 1 

-0 

-0 0 

-G 

0 

1 

-30.0 


.103923 

1 . 

E7 1 . E7 

.16667 

4.266E6 





1773.964 1696.4334 1544.7604 1325.5739 1048.4535 





P 

2 2 

3 

2 0 1 

-0 

-0 0 

-0 

c 

1 

-30.0 


.103923 

1 . 

E7 1 . E7 

.16667 

4.286E6 





810.9104 640 

.1110 

i 462.2983 279.4207 93 

.4816 





P 

3 3 

4 

2 0 1 

-0 

-0 0 

-0 

0 

1 

-30.0 


.103923 

1 . 

E7 1. E7 

.16667 

4.286E6 





-93.4816 -279.4207 -462.2983 - 

■640.1110 -810.9104 





P 

4 4 

5 

2 0 1 

—0 

-0 0 

-0 

0 

1 

-30.0 


.103923 

1 . 

E7 1. E7 

.16667 

4.286E6 






-1048.4535-1325.5739-1544.7604-1696.4334-1773.964 

L 
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8<,6 Example 6: Buckling of a Plate Subjected to Inplane Bending and 

Compression 

The data for this example problem is shown in Figure 8.5. The plate 
is simply supported along all exterior edges. For this problem, M 
is applied initially and then the critical load, corresponding 
to P of the figure, is computed. All elements are divided into 5 
subelements and the N-j-j loads simulating the initial bending moment, 
M, are applied to each subelement as shown. The magnitude of the 
initial bending moment is such that it results in a load of 
46„885 IbSo/in. at the sides of the plate. 
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E = 1 x 10' Ibs/in. 
v = .3 

G = 3.85 x 10 6 Ibs/in. 
t = .05 in. 


2.5" 





Figure 8.6 Plate Subjected to Inplane Moment and Compression 
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PROGRAM S 0 3 4 4 A / BlJ CLASP 3 


NOV 01 72 


THERMAL stress AND BUCKLING Of LAMINATED STIFFENED PANELS 
LOADING -- IN-PLANE BIAXIAL , THERMAL 

BOUNDARY CONDITIONS -- (1) SIMPLY SUPPORTED ALONG EDGES PARALLEL TO THE TRANSVERSE AXIS 

(2) ELASTICALLY RESTRAINED ALONG ANY EDGE PARALLEL TO THE LONGITUDINAL AXIS 
TYPES OF ELEMENTS — (i ) FLAT PLATE 

(2) CURVED PLATE 

(3) BEAM 

EXAMPLE NO. 6 FLATE SUBJECTED TO INPLANE MOMENT AND OPPRESSION 


THIS IS AN EDITTED EXAMPLE OF THE PROGRAM OUTPUT. 
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program control options 


BUCKLING LOA D AND RELATIVE DISPLACEMENTS 

NUMBER OP ITERATIONS ALLOWED IN BUCKLING CALCULATION = 100 

NUMBER or ITERATIONS USED IN EIGENVECTOR CALCULATION = 4 


PROBLEM CHARACTERISTICS 

NUMBER OR NODES = 5 

NUMBER cr FLAT PLATE ELEMENTS = 4 

NUMBER OF CURVED PLATE ELEMENTS - 0 

NUMBER or BEAM ELEMENTS ' 0 

PANEL LENGTH = 15.0GG 

ANALYSIS OPTIONS 

LOAD OPTION (LOAD CP) = ? 

WAVE NUMBER SEARCH CfTlON (MOFT) = 2 

INITIAL VALUE CP LONGITUDINAL WAVE NUMBER (MMl ) = 2 

FINAL VALUE CP LONGITUDINAL WAVE NUMBER (MMA ) = 3 

LOWER LIMIT FOR ROOT SEARCHING CRITERIA (U.W) = 0 

UPPER LIMIT FOR ROOT SEARCH I NG CRITERIA (NUP) = 1 
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NODAL DATA 


- 

coordinates 

BOUNDARY CONDITIONS 


SFRING STIFFNESSES 


NOCE 

Y Z 

CODE 

KCW) 

K (THETA) K (V) 

K (U) 


1 

0 .00000 

0 .00000 

1 

2 

2.50000 

0.00000 

0 

3 

5 .00000 

O.OOGOG 

0 

4 

7.500GQ 

0.00000 

0 

5 

to . 00000 

0.00000 

1 

SUB 

STRUCTURE 

D E F I N I 

T I O N DATA 


mjmbef of substructures =. 1 

substructure: first ncce last node 

START 1 5 
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FLAT PLATE ELEMENT US Eft NUMB Eft 1 USER 

NUMBER OF BUCKLING SUBELEMENTS = 

NUMBER OF RELATIVE DISPLACEMENT SUBDIVISIONS = 
MATERIAL FftOPEftTY INPUT OPTION = 

MATERIAL PROPER TY INPUT FCft MAT OPTION = 


LAYER 

THICKNESS 

EXX 

1 

S.OOOOGOE-02 

1 .0GGGG0E+07 

A- MATH I X 


549450.549 

164835.165 


164835.165 

549450.549 


0 .000 

O.GDO 

B- MATRIX 


-.000 

-.000 


-.000 

-.000 


0.000 

0.000 

D- MATRIX 


114.469 

34.341 


34.341 

114.469 


0.000 

0.000 

INPUT 

LOADS 



NODE I 1 . USER NODE J 2 


EYY NUXT nutx 

.0000DCE+07 3.0000DQE-01 3.00GGQGE-G1 


0 .ODD 
Q.OGQ 
19250G.DGQ 


0 .000 
O.GOQ 
-.000 


Q .000 
O.GGO 
40.104 


SUBSTRUCTURE SUE EL EM ENT 

START 1 

2 

3 

4 

5 


Nil fs£2 

-4.454075E+01 -0. 

-3.985225E+01 -0. 

-3.516375E+G1 *0. 

-3.047525E+01 -0 . 

-2.578675E+01 *0. 


G 

3.85GGGGE+G6 
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PLATS 


ELEMENT 


NET OFFSETS 


ELEM. 


NO- 

1 

2 

3 

4 


START 20 
G.GGGOG 
O.GGGOO 
O.QGGGO 
0 ,GGGGO 


START YO 
G.GGGGO 
G.GOOGO 
O.GGGGO 
O.OGGGO 


END ZO 
G .GGOGO 
O.GGGGG 
O.GGGGG 
O.GGGGG 


END YO 
O.GGGGG 
O.GGGGG 
C.OGGGO 
O.OGOOG 
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E L E 

MEN 

T T R A 

N S F 0 RM A T I 0 

N DAT 

A 



ELEM. 

NO. 

WIDTH 

transformation 
SIN(I) COSU) 

SIN(J) 

COSCJ) 

EL. 

TTFE 


1 

2.50000 

0.00G00 

1 .GGQGO 

0 .00000 

1 .00000 

FLAT 

PLATE 


2 

2 . 50000 

0.00000 

1 .00000 

0.00000 

1 .00000 

FLAT 

FLATE 


3 

2.50000 

0.00000 

1.00000 

0.00000 

1 .00000 

FLAT 

PLATE 


4 

2.50GG0 

0.00000 

1 .00000 

0.00000 

1 .GGGGG 

FLAT 

FLATE 


ELEMENT SUBSTRUCTURE ASSIGNMENT CATA 

NUMBER BLOCK 

1 START 

2 START 

3 START 

4 START 
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4 * 

* MOCC IS M = 2 * 

* * 

* * 


*4** *«$ $4$ v 


EXAMPLE NO. 6 PLATE SUBJECTED TO INPLANC MOMENT AND COMPRESSION 


VARIABLE LOAD CM 

number cf neg. 

EUCKLING DETERMINANT 

(A* 

FIRST PLATE ELEMENT 

CM DIAGONAL 

A 

6 

50.00000000 

1 

-.748 72 

148 

25.00000000 

0 

.46735 

152 

ST. 50000000 

0 

.18799 

152 

43.75000000 

0 

.06517 

152 

46.87500000 

0 

.12550 

148 

48.43750000 

1 

-.31694 

148 

47.65625000 

1 

-.09706 

148 

47.26562500 

0 

.22216 

144 

47.46093750 

1 

-.66673 

144 

47.36328125 

1 

-.22262 

144 

47.31445313 

1 

-.08064 

136 

47.29003906 

c 

.11090 

144 

47.30224609 

0 

.88460 

140 

47.30834961 

0 

.43976 

140 

47.31140137 

0 

.21736 

140 

47.31292725 

0 

.10616 

140 

47.31369019 

0 

.80893 

136 

47.31407166 

0 

.36414 

136 

47.31426239 

0 

.14175 

136 

47.31435776 

0 

.48894 

132 

47.31440544 

1 

-.40062 

132 

47.31438160 

0 

.70657 

128 

47.31439352 

1 

-.17823 

132 

47.31438756 

1 

-.06703 

132 

47.31438458 

1 

-.18298 

128 

47.31438309 

0 

.26179 

128 

47.31438383 

0 

.63041 

124 

47.31438421 

1 

-.07179 

128 

47.31438402 

1 

-.25915 

124 


COMMENTS 
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CLCHCMT load s 


tLt>, 


fO. 

Nil 4 

N22 

EFSU 

EFS22 

1 

47.314384 

0 . 

9.462876806-05 

-2.636863046-05 

z 

47.314384 

0 . 

9 .4628 7680 £-05 

-2.63886304E-05 

3 

47.314384 

0 . 

9 .<628 7680 £-05 

-2 .838863046-05 

4 

47.314384 

0 . 

9 .<628 7680 £-05 

-2 .83686304E-05 

TOTAL 

AXIAL 

TR A*I5 VU\S£ 

load 

AT DUCKLING 
473 

0 . 




P~D£Av 
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EXAMPLE NO. 6 PLAT 

euckltno lo 

CRITICAL WAVE NUMEER M 

H 

Z 


r SUCJECTEC TO I NrLAflC MOMENT AfC CCM^ESSlCN 
AC COMPUTATION summart 

= 2 
LCAD 

47 . 3143840 ? 




EXAMPLC NO.6 PLATE SUBJECTED TO INPLAHC MOMENT AND COMPRESSION 


i * 4 9 « ♦ * * 4 4 * «> * $ fc* £ £ * * *4 ** *«> •> *> $ * 4 


EIGENVECTOR ♦ 

AND A 

RELATIVE DISPLACEMENTS * 

***<t*<U***£**l*J>l*£«* ££0 3 £*$£«**££« 

M = 2 4 


l £*£*££££££££££££££££££££££«£££ 


STIFFNESS MATRIX BANDWIDTH = 6 


EIGENVECTOR 


ITERATION NUMBER 
normalizing FACTCR 


1 

3.4 7454 71 E + GO 



SUBSTRUCTURE 

NODE 

W 

theta 

V 

START 

2 

-8 . 99468 629E-0X 

3. 036841 42E-02 

4.20628366E-01 


3 

-1 .OOGGOOOGE+OO 

3.63699823E-01 

4.179 5234 7E-01 

EIGENVECTOR 

4 

-5.01184145E-01 

2.87807295E-01 

2.878G 7438E-01 


U 

4.00709375E-01 
3 .348418 56£-Ol 
1 . 209 62283 E~06 


ITERATION number 2 

NCRMALIZING FACTOR 4 .3771 381 E+07 


SU&STRUCTCRE NODE 


W 


THETA 


V 


u 


START 2 

3 

4 


EIGENVECTOR 


8.13892665E-01 
1 .CGOCCGQOZ+OO 
6.0794 75 59E-G1 


2.19988135E-01 
-6.389306S6E-02 
”2.21 756339 £-01 


7. 4 63 54332 £~ 14 
9.308641 67E- 14 
8.67101555E-14 


4.83778C54E-14 
3.78 62Q706E-14 
1 .05354926E-1 5 


ITERATION NUMBER 3 

NORMALIZING FACT OR 4 .41 1 31 57£*07 


SUeSTRUCTLRE NCCE 


W 


THETA 


V 


U 


START 2 

9 
4 


EIGENVECTOR 


-0.13892665£-Ol 
-I .GGGOGGGG£*GO 
-6.07947559E-0I 


-2.19968135E-0I 
6.3893G667E-02 
2 .21 766339£-0l 


1 .62043499E-26 
2.03044411 £-26 
1 .99944678E-26 


G . 56356820^-2 7 
4 .91129494E-27 
-4.01879091 E-30 


ITERATION NUMCCR 4 

NCR HA L IZ I NG FACTOR 4 .4 1 1 3 1 5 7£* 0 7 

SUSS TR UC T tR E NOCC W THCTA 


8.44 



START 


4 

z 

3 


•.136926656-Gl 
1 .OOOOOOCQr+OO 


2 .199881356-01 
-6.389306S7E-02 
-2 .2! 766329 r_ 0! 


3.680903166-39 

4.447875756-39 


9 .99 70 5 5636-40 
6.878105036-40 
-8 . '232?! 70 r - ii 


8.45 



.EXAMPLE no. 6 PLATE SUD JECTCC TO lNPLAHC MOMENT AM CCHF^CSSICN 


4 4 *44 *>*44 44 4 4444 4 44 * 


4 4 

4 RELATIVE CISPLACEMENTS 4 

* A 
<»**;>:*<ttt4:i:*:*44444444444444*44444 

* * 

4 M = 2 4 

4 4 

4«.^44444444444444^44444444 44444 


START SECTION — BLCCK W>. 1 


ELEMENT NO.= 

1 

TYPE = 

FLAT 

plate nice 

I 

= 1 NICE 

J - 

2 

WICTH = 

2.5000 

r 


w 


V 


theta 





-.2500 


.GOG DO 


.00000 


.38339 





-.1250 


.23715 


.00000 


.37160 





0.0000 


.45967 


.DOGGO 


.33768 





.1250 


.65542 


.00000 


.23 541 





.2500 


.81369 


.00000 


.21999 





ELEMENT NI>.- 

2 

TYFE = 

FEAT 

PLATE NICE 

I 

= 2 NICE 

J = 

3 

WICTH = 

2.5000 

r 


W 


V 


THETA 





-.2500 


.81389 


.00000 


.21999 





-.1250 


.92866 


.00000 


.14716 





0.0000 


.99747 


.00000 


.07250 





.1250 


1 .02013 


.00000 


.00058 





.2500 


1.00000 


.00000 


-.06389 





ELEMENT N9- = 

3 

TYPE = 

FLAT 

PLATE NjCE 

I 

= 3 NICE 

J = 

4 

WICTH = 

2.5600 

T 


W 


V 


THETA 





-.2500 


1 .0 0000 


.00000 


-.06389 



t 


-.1250 


.94223 


.00000 


-.11926 





0.0000 


.85316 


.00000 


-.1.6394 





.1250 


.73955 


.00000 


-.19783 





.2500 


.60795 


.00000 


-.22177 





ELEMENT NO. = 

4 

TYPE = 

flat 

PLATE NICE 

I 

= < NICE 

J - 

5 

WICTH = 

2.5600 


r 

w 

V 

THETA 

-.2500 

-60795 

.00000 

-.22177 

-.1250 

.46410 

.00000 

-.23728 

0.0000 

.31270 

.00000 

-.24 629 

.1250 

.15720 

.00000 

-.25066 

.2500 

.00000 

.00000 

-.25194 


8.46 



9.0 PROGRAM DESCRIPTION 

This section is a description of the organization and function of the 
various routines included in BUCLASP3. 


9.1 Overlay Structure 

BUCLASP3 consists of a (0,0) level overlay which acts as a monitor in 
selection of primary overlays. 

The overlay structure is: 

(1) (0,0) Overlay 

BUCLASP ,0 ,0 S0344A, Monitor 


(2) Primary Overlays 

(a) BUCLASP, 1,0 DATAPRO, Reads all user input and outputs 

basic job description. 

(b) BUCLASP, 2,0 THERMAL, Controls calculation of thermal 

displacements and stresses 

(c) BUCLASP, 3,0 LOADING, Generates element stiffness matrices, 

merges panel stiffness matrix, and controls 
calculation of panel buckling load. 


(d) BUCLASP, 4,0 DISPLAC, Calculates for the panel buckling 

load, the panel eigenvector and element 
relative displacements. 


(3) Secondary Overlays 

(a) BUCLASP, 2,1 TSTIFF, Generates panel thermal stiffness and 

loads matrices 

(b) BUCLASP, 2, 2 TDISPL, Merges loads and calculates global 

thermal displacements. 


9.1 



(c) BUCLASP.2 ,3 


TSTRESS, Computes element stresses and 
complementary displacements. 


9.2 



Main Overlay 0,0 
Program Name: S0344A 

Subroutines: ELMDEF, ELMFCH, ERROR, IGAL, KELNUM, PAC, PRINT, UNPAC 


Primary Overlay 1.0 Primary Overlay 3.0 
Main Program: DATAPRO Main Program: LOADING 


Subroutines: RDVLS 

USERID 
BLKDEF 
PLTDEF 
TBPOINT 
NEXTC 
FINDIN 
CLNODE 
RDTBLE 
IBLKNO 
STRMOV 
TDOUT 
SDOUT 
BLKELM 
SORT 
TDATA 


Subroutines: BLKWRT 

LCNTRL 
DM 

PLATE 

SBKDT1 

UPPRBD 

UPRBND 

DE 

DH 

SPGCNS 

STRFR 

SYMDET 

GAPUPP 

DB 

BLKGEN 

TRANF 

SBLKDT 

SELIM 

REDUVE 

DBLERT 

STORE 

PRDMTX 

PROOT 

HP LATE 

MOVX 

MLTPLY 

FPLATE 

FPLCOL 

BEAM, 

SOLVV 

DT 

RGEN 

CDTM 

VIPDR 

ZARK 

MATZ 

TRANSF 


Primary Overlay 4.0 


Main Program: DISPLAC 


Subroutines: BANDW 

COMPAC 
BLKRED 
SYMBND 
REDUVE 
DIS 

FETCHO 

READTR 

FNLDIS 

TRAND 
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Secondary Overlay 2.1 ■ 

Secondary Overlay 2.2 

Secondary Overlay 2.3 

Main Program: TSTIFF 

Main Program: TDISPL 

Main Program: TSTRESS 

Subroutines: TRAND 

Subroutines: TSTORE 

Subroutines: TRAND 

TPLATE 

BANDW 

DTMOD 

YY 

COMPAC 

LI NIC 

YB 

SYMBND 

FTCHDF 

TBEAM 

REDUVE 



SMOD 

TRHS 

FMODB 

SPRING 

PLACE 

TRANB 

SPROOT 

DT 

SRGEN 

CDTM 

VIPDR 

SPLATE 

SBEAM 

SPLCOL 

SOLVE 

ZARK 








9.2 Basic Program Flow 



9.5 












DATAPRO 



9.6 













yes 







Return 




9.8 







9.9 











9.3 Program/Subroutine Description 
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lt*6 


Routine 

Purpose 

Configurations 

Tests 

BANDW 

BANDW calculates the bandwidth of 

Panel with 

F-l 


the panel stiffness matrix for 

1) 

Start block only 

F-38 


COMPAC. 

2) 

Start and end blocks only 

F-8 



3) 

Start, Repeat and End blocks 

F-22 

BEAM 

BEAM generates stiffness matrices. 

1) 

Beam element is input at an angle 

F-23 



2) 

Beam element is input with offsets 

F-39 

BEAMDF 

BEAMDF reads and processes beam 

1) 

General beam element 

F-25 


material properties 

2) 

Rectangular beam with 





a) One lamina 

F-24 




b) More than one lamina 

F-23 



3) 

Circular beam with 





a) One lamina 

F-39 




b) More than one lamina 

F-26 



4) 

Laminated beam input by 





a) B-i cards 

F-23 




b) Tables 

F-39 
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Routine Purpose Configurations Tests 

BLKDEF BLKDEF processes block definition node 1. Problem with 

pairs to produce the row and column a) with equivalent nodes F-23 

definition for each block and to b) without equivalent nodes F-22 

produce the corresponding merging c) start block only F-l 

information. d) start and end blocks only F-38 

e) start, repeat and end blocks F-22 

f) beam elements F-23 

BLKELM BLKELM constructs a list of all elements Problem with 

in a substructure. 1. start, repeat and end blocks F-12 

2. plate element F-l 3 

3. beam element F-23 

BLKGEN BLKGEN forms the buckling determinant 1. Problem with F-l 1 

for a given load and wave number m. a) varying F-l 5 

b) N-^ with biaxial ratio F-28 

c) N 22 varying F-l 9 

d) N 22 varying with constant strain F-20 

2. Problem with 

a) Clamped Node F-2 

b) Simply supported Node F-l 

c) Free Node F-l 

d) Sprung Node F-l 1 

3. Problem with plates and beams F-23 



er« 


Routine 

Purpose 

Configurations 

Tests 



5o LOADOP = 5 

B-l 



6. LOADOP = 6 

B-2 



7. LOADOP = 7 

B-3 



8. LOADOP = 8 

B-4 

BLKRED 

BLKRED reads from disk file repeat 
substructure of panel matrix 


F-22 

BLKWRT 

BLKWRT writes repeat substructure of 
panel matrix to disk file. 


F-22 

CDTM 

CDTM evaluates the determinant of 




a complex square matrix 


F-l 

CLNODE 

For each node between the first 

1. Plates only 

F-l 


node and the last node of a block 
definition pair, investigate the 
elements having the node as a 
definition node. From those elements 
find the smallest and largest 
intern or node number. 

2. Plates and Beams 

F— 23 



frt*6 


Routine 

Purpose 

Configurations 

Tests 

COMPAC 

COMPAC takes the panel stiffness matrix 

Panel with 



as generated in the buckling computa- 

a) Start block only 

F-l 


tion and puts it into symmetric band 

b) Start and end blocks only 

F-38 


format suitable for eigenvector 

c) Start, Repeat and end blocks 

F— 22 


computation with all repeat blocks 

d) Start, Repeat and end blocks in 

F-10 


present. 

curved panel option 


DATAPRO 

DATAPRO read input cards 1-9, cards 

Problem with 



T, C, E, D, N. S, P, B, F, and L and 

1 . Start block only 

F-l 


process the corresponding data. 

2. Start and end blocks only 

F-38 


The remaining data processing routines 

3. Start, Repeat and End blocks 

F-22 


are called as required. 

Problem with presence of 




4. T cards 

F-39 



5. C cards 

F-23 



6. E cards 

F-23 



7. D cards 

F-23 



8. N cards 

F-23 



9. S cards 

F-23 



10. B cards 

F-23 



11. P cards 

F-23 



12. F cards 

F-22 



••15 


Routine 

Purpose 

Configurations 

Tests 



13. L cards 

F-22 



Problem with 




14. Beam element 

F-23 



15. Flat plate element 

F-23 



16. Curved plate element 

F-21 



User Diagnostics 

F-40 - F-90 

DB 

DB controls the evaluation of the panel 
buckling load 


F-l 

DBLERT 

DBLERT is the error trap routine for 
PROOT encountering a double root. 


Inspection 

DE 

DE controls merging and buckling 

Problem with 



condition determination of 

1. LOADOP = 6 or 8 

B-2 


subelement approximation of an 

2. LOADOP = 5 or 7 

B-3 


element. 

3. 2 subelements per element 

B-l 



4. 3 subelements per element 

B-9 



5. Offsets at left end of element 

5-5 



6. Offsets at right end of element 

5-5 



91*6 


Routine 

Purpose 

Configurations 

Tests 

DH 

DH generates and evaluates the 

Thermal problem with 



determinant of the retained nodal 

1 . 

L0AD0P = 5, 7 

B-3 


stiffness of a plate. 

2. 

LOADOP = 6, 8 

B-4 



3. 

Plate with 





a) Simply supported node 

B-9 




b) Clamped node 

Inspection 




c) Free node 

B-9 




d) Sprung node 

Inspection 

DIS 

DIS computes the relative displacements 

1 . 

Panel with 



of the deformed panel . 


a) Start block only 

F-l 




b) Start and End blocks only 

F-38 




c) Start, Repeat and End blocks 

F-22 



2. 

Problem with 





a) Plate with zero B-matrix 

F-29 




b) Plate with non-zero B-matrix 

F-7 

DISPLAC 

DISPLAC controls the execution BUCLASP, 



Inspection 


3,0 overlay which computes panel 



F-l 


eigenvector and element relative 
displacements. 
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Routine 


Purpose 

DM DM generates and evaluates the 

determinant of the retained node 
on a plate with a non-interior 
node. 


DT DT calculates the determinant of 

the coefficient matrix for the 
equilibrium equations. 

DTMOD DTMOD performs backsubstitution for 

deflection on reduced out nodes. 


Configurations 


Tests 


1. a) 

N i i varying 

NF-11 

b) 

Nil varying with non-trivial 



biaxial ratio 

F-l 5 

c) 

N 22 varying 

F-28 

2. Plate with 


a) 

Simply supported node 

F-2 

b) 

Clamped node 

F-l 

c) 

free node 

F-l 

d) 

Sprung node 

F-l 1 

1. FI ate plate with non-zero B-matrix 

F-7 

2. FI ate plate with zero B-matrix 

F-29 


Thermal plate problem with boundary 
condition 


S-l 



8t*6 


Routine 

Purpose 

Configurations 

Tests 

ELMDEF 

ELMDEF retrieves for an element in a 

Problem with 



substructure its type, global number. 

1 . 

Plate in 



local number and plate lamina number. 


a. Start Substructure 

F-22 




b. Repeat Substructure 

F-22 




c. End Substructure 

F-22 



2. 

Beam in 





a. Start Substructure 

F-22 




b. Repeat Substructure 

F-22 




c. End Substructure 

F-22 

ELMFCH 

ELMFCH retrieves from element list 

Problem with 



element type, plate lamina number 

1 . 

Beam element 

F-22 


and element local number given 
element global number. 

2. 

Plate element 

F-22 

ERROR 

ERROR prints name of routine where 



Inspection 


error was detected and error number 



F-40 


and calls exit. 




EVALUE 

The routine EVALUE extracts the 
lowest eigenvalue from a general 
eigenproblem of up to order 15 for 
the GALERKIN method. 



F-17 



6r6 


Routine 

FETCHD 


FINDIN 


Purpose 


Configurations 


Tests 


From the panel eigenvector, fetch 
the nodal deflection component 
for an interior node. 


a) For the first node in a block 
definition pair, find the 
smallest interior internal node 
larger than the first node. 

b) For the last node in a block 
definition pair find the 
largest interior internal node 
smaller than the last node. 


1 . Node is in row 
a) Start block 

a.l) Is part of an element with a F-22 
node in the repeat block. 
(Interrelationship elements 
defined) 



a. 2 Is not 

F-22 

b) 

Repeat block 

F-22 

c) 

End block 

F-22 



F-22 
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Routine 

Purpose 

Configurations 

Tests 

FMODB 

FMODB and FMOD produce final plate 

Problem with thermally loaded plate with 



thermal loads for merging 

1 . 

left node constrained 

S-l 



2. 

Right node constrained 

S-l 



3. 

Both nodes unconstrained 

S-l 

FNLDIS 

For a beam or plate element obtain 

1 . 

Beam element 

F-23 


the element displacement vector by 

2. 

Plate element 



fetching the element nodal 


a) left node reduced out 

F-28 


displacement vector or backsubsti- 


b) right node reduced out 

F-28 


tuting for reduced out nodal 
components of plate elements. 

The W-j's are solved for plate elements. 


c) both nodes interior 

F-l 

FPLATE 

FPLATE (1) calculates the P-root 

1 . 

Plate element with non-zero B matrix 

F-l 


normalization factor for the set 

2. 

Plate element with zero B matrix 

F-29 


of P roots for a plate element, 
(2) puts together the deflection 
and force matrices for the plate. 
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Routine 


Purpose 

FPLCOL FPLCOL generates for one P-value the 

columns of the deflection and force 
matrices for the left and right side 
of a plate element. 


Configurations Tests 

1. Curved plate element with non zero 

B-matrix. varying F-31 

2. Curved plate element with zero 

B-matrix. varying F-30 

3. Flate plate element with non zero 

B-matrix. N-. -j varying F-32 

4. Flate plate element with zero 

B-matrix. varying F-34 

5. Curved plate element with non zero 

B-matrix. N«£ varying F-7 

6. Curved plate element with zero 

B-matrix. varying F-37 

7. Flate plate element with non zero 

B-matrix. varying F-29 

8. Flate plate element with zero 

B-matrix. varying F-28 
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Routine 

Purpose 

Confi gurations 

Tests 

FTCHDF 

FTCHDF retrieves from the panel 

Thermal problem with node in 



thermal displacement vector the 

1 . 

Start substructure 

S-4 


nodal deflection component for an 

2. 

Repeat substructure 



interior node. 


a. Element connecting to start 





substructure 

S-4 




b. Element connecting to end 





substructure 

S-4 



3. 

End substructure 





a. Element connecting to repeat 





substructure 

S-4 




b. Element not connecting to 





repeat substructure 

S-4 

GAPUPP 

GAPUPP uses Galerkin's Method to 

1 . 

Flate plate 



obtain an approximation to an 


a) N-|-| varying 

F-l 1 


individual plate's critical load. 


b) N 22 varying 

F-28 




c) Biaxial ratio on 





1) first plate element 

F-l 5 




2) subsequent plate element 

F-16 



2. 

Curved plate 

F-16 



cr 6 


Routine 

Purpose 

Confi gurations 

Tests 

HPLATE 

HPLATE produces the deflection. 

Problem with 



force and stiffness matrices for 

1 . 

Curved plate 

F-21 


a subdivided plate element. 

2. 

Flat plate 

F-22 



3. 

Offsets at left end 

S-5 



4. 

Offsets at right end 

S-5 



5. 

LOADOP = 5 

F-ll 



6. 

LOADOP =2,3 

F-28 



7. 

LOADOP =5,7 

B-l 



8. 

LOADOP =6,8 

B-3 

IBLKNO 

Obtain for an interior, internal node 

1 . 

Start block only 

F-l 


the corresponding row block number. 

2. 

Start and End block only 

F-38 



3. 

Start, Repeat and End blocks 

F-8 

I GAL 

IGAL supplies information as to the 
length of blank common in each overlay. 



Inspection 
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Routine 

Purpose 

Configurations 


Tests 

KELNUM 

KELNUM calculates for an element its 

Thermal problem with 




local and global thermal element 

1. 

Beam in 




number. 


a. Start structure 


S-5 




b. Repeat structure 


S-4 




c. End structure 


S-5 



2. 

Flat plate in 






a. Start structure 


S-4 




b. Repeat structure 


S-4 




c. End structure 


S-4 



3. 

Curved plate in 






a. Start structure 


S-5 




b. Repeat structure 


S-5 




c. End structure 


S-5 

LCNTRL 

LCNTRL finds for an interval of real 

1. 

No buckling load found 


Inspection 


numbers the buckling load (if it 

2. 

Buckling load found 




exists) in the internal to a specified 

3. 

Double root encountered at a 

trial 

load 


tolerance using a bisection technique. 

4. 

Zero determinant encountered 

at a 





trial load. 
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Routine 


Purpose 

LININC L IN INC controls print line for thermal 

stress output. 

LOADING LOADING controls the computation 

of the critical load and wave 
number for the panel . 


MATZ MATZ zeros out a specified block of 

rows and columns of a matrix. 

MLTPLY MLTPLY performs product of two 8X8 

matrices. 

MOVX MOVX transfers force and deflection 

matrices. 


Confi gurations 


Tests 


S-l 


1. Various wave number search options 

a) 1 F— 37 

b) 2 F-l 

c) 3 F-36 

d) 4 F-39 

2. Upper Bound Input F-37 

3. Upper Bound Calculated. F-l 1 

4. Upper Bound Calculation Option Only F-l 1 

5. Problem with Beams, FI ate and Curved F-40 

Plates 


F-l 


B-l 


B-l 
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Routine 


Purpose 

NEXTC NEXTC returns next character on 

table pointer card. 

PAC will pack the (J-I+l) 
rightmost bits of the word B 
into the bit positions I-J 
of the word A. The remaining 
bits of A are unchanged. 

PLACE PLACE merges a 4 X 4 element 

submatrix associated with a node 
pair into the appropriate sub- 
block of the panel thermal 
stiffness matrix. 

PLTDEF PLTDEF reads and processes plate 

definition data to produce plate 
lamina extensional, coupling and 
bending stiffness matrices. 
PLTDEF produces transformation 
matrices used in strain and 
load calculation. 


PAC (A, I, 
J, B) 


Configuration 

1. Test for end-of-file 

2. Illegal character 


Tests 


F— 37 
F-l 


Thermal problem with 

1. Start substructure only S-l 

2. Start and End substructures only S-6 

3. Start, Repeat and End substructures S'-5 

Problem with 


Lamina stress-strain input 


a) by P-2 card 

F— 37 

b) by Table 

F— 37 

Engineering constants input 


a) by P-1 card 

F-l 

b) by Table 

F— 37 

Non-zero B matrix 

F-30 

a) have same radius, loads and 

F-40 

sti ffness 
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Routine 

Purpose 

Configuration 

b) have different radius, loads and 
stiffness 
4. Zero B matrix 

Tests 

F-32 

F-33 

PRDMTX 

PRDMTX reduces out the effects of a 
constrained node for a plate element. 

1. The element nodal intercoupling 
matrix is singular 

2. The element nodal intercoupling 
matrix is non singular 

Inspection 

F-l 

PRINT 

PRINT prints a rectangular matrix up 
to 8 columns at a time and all rows. 


F-81 

PROOT 

PROOT calculates the roots of the 
equilibrium equations for a plate 
element. 

1. Plate element with non-zero B-matrix 

2. Plate element with zero B-matrix 

F— 27 
F— 27 

RDTBLE 

Reading from card images (columns 
11-80) floating point data in a 
free-field format with commas and 
end of physical card as field 
delimeters. Slash (/) is the 
logical card delimeter limit of 
fourteen significant digits. 

1. 2 or more physical cards used 

2. E format data 

3. E format data with 14 signicant digits 

4. F format data 

5. F format data with 14 signicant digits 

6. Error exits 

a) end-of-file, missing slash 

F— 37 and 

subroutine 

tests 
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Routine 


Purpose 


RDVLS Reads pairs of integers in a 

2(2X, 13) format from card images 
and stores non zero pairs in an 
array. 

READTR Read from a binary file a nodal or 

elemental transformation matrix. 

REDUVE Reduces one vector into another 

vector 

COMPASS Routine 

RGEN RGEN generates the elements of the 

coefficient matrix for the 
equilibrium equations for the 
current plate element. 


Configuration 


Tests 


b) Illegal data, bad character, too 
many characters. 

c) Table size limit exceeded 

d) No data input 

1. Eof of file encountered F-l 

2. Maximum of receiving array is 
exceeded. 

1. Nodal transformation requested F-l 

2. Elemental transformation requested F-l 

F-l 


1. Curved plate element with zero 

B-matrix. N,-j varying F-7 

2. Flate plate with non-zero B-matrix. 

N-|i varying F-37 

3. Flate plate with zero B-matrix. 

N-| i varying 


F-29 



9.29 


Routine 


Purpose 


SBEAM SBEAM generated beam element 

stiffness matrix 


SBLKDT SBLKDT computes the determinant of 

a symmetric matrix which has an 
overlapping and repeated block 
structure (where the blocks differ 
in content) 


Configuration Tests 

4. Curved plate with non-zero B-matrix. 

N-|-| varying F-28 

5. Curved plate element with zero 

B-matrix. varying F-31 

6. FI ate plate with non-zero B-matrix. 

N 22 varying F-30 

7. FI ate plate with zero B-matrix. 

^2 varying F-32 

8. Curved plate with non-zero B-matrix. 

^2 varying F-34 


1. Thermal beam element is input at an 


angle S-5 

2. Thermal beam element is input with 

offsets Inspection 

Symmetric Matrix with 

1. Start block only F-l 

2. Start and end blocks only F-38 

3. Start, 1 repeat and end blocks F-8 

4. Start, 2 repeat and end blocks F-22 



9.30 


Routine 

Purpose 

Configuration 

Tests 

SBKDT1 

SBKDT1 calculates the determinant of 
a symmetrix matrix for upper bound 
calculation. 



F-17 

SDOUT 

SDOUT outputs to scratch disk the 

Externally loaded problem with 



load data by substructure. 

1 . 

Start substructure only 

B-l 



2. 

Start and end substructures only 

B— 22 



3. 

Start, repeat and end substructures 

B— 21 



4. 

Plate element 

B-21 



5. 

Beam element 

B-21 

SELIM 

SELIM performs a Gaussian reduction 

1 . 

Non-zero pivot 

F-l 


on a symmetric matrix. 

2. 

Zero pivot 

Inspection 

SMOD 

SMOD performs steps of reduction of 

Thermal problem with plate with a 

S-l 


plate elemental stiffness with 

boundary condition 



constrained node. 




SOLVE 

SOLVE solves the matrix problem 

1 . 

X T is non-singular 

S-l 


X-A = B for X 

2. 

X is singular 

Inspection 



lf*6 


Routine 

Purpose 

Confi guration 

Tests 

SOLVV 

Solves matrix problem X.A = B for X 

1 . 

X T is non-singular 

F-l 


or A«X = B. 

2. 

X^ is singular 

Inspection 



3. 

X is singular 

B-l 



4. 

X is not singular 

Inspection 

SORT 

SORT sorts list of integers from 
low to high. 



B-l 

SPGCNS 

SPGCNS produces the spring constants 

Problem with LOADOP = 5, 6, 7, 8 where 



for the end subelements of an element. 

element has 




1 . 

interior node 

B-9 



2. 

simply supported node 

B-9 



3. 

clamped node 

Inspection 



4. 

free node 

B-9 



5. 

sprung node 

Inspection 

SPLATE 

SPLATE calculates the P-root 

1 . 

Plate element with non-zero 



normalization factor for the set 


B-matri x 

S-l 


of P roots for a plate element and 
puts together the deflection and 
force matrices for the plate. 

2. 

Plate element with zero B-matri x 

S-2 



9.32 


Routine 


Purpose 

SPLCOL SPLCOL generates for one P-value 

the columns of the deflection and 
force matrices for the left and 
right sides of a plate element. 


SPRING SPRING generates the nodal spring 

stiffnesses for a plate element. 


SPROOT SPROOT calculates the roots of the 

equilibrium equation for a plate 
element in a thermal environment. 


Tests 


Conf i guration 
Thermal problems with 


1. Curved plate element with non-zero 
B-matri x 

S-2 

2. Curved plate element with zero 
B-matri x 

S-3 

3. Flat plate element with non-zero 
B-matrix 

S-7 

4. Flat plate element with zero 
B-matrix 

S-l 

Thermal problem with plate with 
1. Both nodes interior 

S-l 

2. Free node 

Inspection 

3. Clamped node 

Inspection 

4. Simply supported node 

S-l 

5. Sprung node 

S-3 


1. Plate element with non-zero B-matrix S-l 

2. Plate element with zero B-matrix S-2 

3. Plate element which is repeated. S-l 



9.33 


Routine 


Purpose 

SRGEN SRGEN generates the elements of the 

coefficient matrix for the 
equilibrium equations for the 
current thermal plate element. 


STORE STORE merges a 4 X 4 element submatrix 

associated with a node pair into the 
appropriate subblock of the panel 
stiffness matrix. 

STRFR STRFR stores nodal stiffness 

components into compact storage 
scheme for element. 

STRMOV Moves characters starting at a 

position in one array to a different 
position in another array. 


Configuration Tests 

1. Curved plate element with zero 

B-matrix S-3 

2. Curved plate element with non-zero 

B-matrix S-2 

3. Flat plate element with zero 

B-matrix S-l 

4. Flat plate element with non-zero 

B-matrix S-7 

Panel with 

1 . Start block only F-l 

2. Start and end blocks only F-38 

3. Start, repeat and end blocks F-22 

Problem with LOADOP = 5, 6, 7, or 8. B-9 


1. Starting array is more than one word Inspection 

2. Receiving array is more than one word F-37 



9.34 


Routine 

Purpose 

Configuration 

Tests 

SYMBND 

Linear equation solver for a symmetric 

1 . 

Zero pivot encountered 

Inspection 


band matrix 

2. 

All non zero pivots 

F-l 

SYMDET 

SYMDET counts the number of neqative 
occuring on the diagonal of a banded 
matrix during decomposition. 



B-9 

S0344A 

S0344A controls the execution of 

1 . 

Multiple Data Sets 



BUCLASP3 through its primary overlays. 

2. 

Data Check Option Only 

F— 41 


provides central processor time 

3. 

Buckling Option only 

F— 36 


usage information and controls 

4. 

Buckling and Eigenvector option 

F-l 


various n/n+1 root crossing 

5. 

Multiple root searching option 

F-35 


searches. 

6. 

Thermal Stress 




7. 

Thermal Buckling 


TBEAM 

TBEAM generates thermal loads for 

Problem with thermally loaded beam 

S-5 


a beam element. 




TBPOINT 

TBPOINT reads the table pointer card 

1 . 

Repeat option left of slash 

F-37 


following P and B cards and constructs 
an array of pointers. 

2. 

Repeat option right of slash 



r* 

k) 



9.35 


Routine 


Purpose 

TDATA TDATA reads and preprocesses thermal 

data for beam and plate elements. 


TDISPL TDISPL controls the merqing of the 

panel thermal right hand side and 
coefficient matrices and the 
computation of the global thermal 
displacements 

TDOUT TDOUT outputs to scratch disk the 

thermal data by substructure 


Confiquration 


Tests 


Problem with 

1 . General beam S-5 

2. Pxectangular beam Inspection 

3. Circular beam Inspection 

4. Plate with 

a. Thermal data in z-direction 

at mid plane S-7 

b. Thermal data in z-direction 

by coefficients S-l 

Thermal problem with 

1. Start substructure only S-l 

2. Start and end substructures only S-6 

3. Start, repeat and end substructures S-5 

4. Plate element S-5 

5. Beam element S-5 

Thermal problem with 

1. Start substructure only S-l 

2. Start and end substructures only S-6 

3. Start, repeat and end substructures S-5 

4. Plate element S-5 

5. Beam element S-5 



9.36 


Routine 

Purpose 

Configuration 

Tests 

THERMAL 

THERMAL controls execution of thermal 

1 . 

Check print option ON 

Inspection 


stress overlays. 

2. 

Check print option OFF 

S-5 

TP LATE 

TPLATE generates thermal loads for a 

Thermal problem with 



plate element. 

1 . 

Flat plate 

S-l 



2. 

Curved plate 

S-2 

TRANB 

TRANB transforms a nodal stiffness 
matrix from the local axis to the 
global axis. 



S-l 

TRAND 

Transforms nodal deflection component 



F-l 


from the global axis to the local 



F-2 


axis. 



F-3 


F-4 

F-5 

F-6 

F-8 

F-9 



9.37 


Routine Purpose Configuration Tests 

TRANF TRANF transforms nodal stiffness F-l 

matrix from the local axis to the F-2 

global axis for plate elements F-3 

with non-interior node. F-4 


F-5 

F-6 

F-8 

F-9 

TRANSF TRANSF transfers specified blocks of F-l 

rows and columns from one array to 
another array. 

TRHS TRHS produces the backsubstitution Thermal problem with plate with a 

matrix for a plate element with boundary condition with thermal loads. S-l 

thermal loads and a constrained 

node. 

TSTIFF TSTIFF controls the generation of the Thermal problem with 


stiffness and 

loads matrix for the 

1 . 

Start substructure only 

S-l 

thermal global 

displacements and 

2. 

Start and end substructures only 

S-6 

the merging of 

the stiffness matrix. 

3. 

Start, repeat and end substructures 

S-5 



4. 

Plate element 

S-5 



5. 

Beam element 

S-5 



9.38 


Routine 

Purpose 

Configuration 

Tests 

TSTORE 

TSTORE merges nodal components of 

Thermal problem. 

S-l 


element right hand side matrices 
into the qlobal right hand side 
matrix. 




TSTRESS 

TSTRESS computes the elemental stresses 

Thermal problems with 



and complimentary displacements. 

1. 

Start substructure only 

S-l 


TSTRESS either (1) calculates 

2. 

Start and end substructures only 

S-6 


thermal stresses at specified points 

3. 

Start, repeat and end substructures 

S-5 


for each element or (2) calculates 

4. 

Plate element 



average and N «2 "for each 


a. Thermal stress mode 

S-l 


subdivision of a plate element or 


b. Thermal buckling mode 

S-5 


the average for a beam element 


c. with zero B-matrix 

S-l 


due to thermal load. 


d. with non-zero B-matrix 

S-2 



5. 

Beam element 





a. Thermal stress mode 

S-l 




b. Thermal buckling mode 

S-5 

IINP AC (A, 

UNPAC extracts bits I-J from the word 



Inspection 

I, J, B) 

A and places the right-adjusted into 
the word B. The remaining bits of 
B are set to zero. 






9.39 


Routine Purpose 

UPPRBD UPPRBD calculates the upper bound 

for the critical load. 


Tests 


Confi guration 

1 . Problems wi th 

a) beam elements F-27 

b) similar plate with 

1) same boundary condition - 

with same and different widths F-33/F-27 

2) different boundary condition F-27 

3) same spring stiffness - 

with same and different widths F-33/F38 

4) different spring stiffness F-38 

c) non neglected plate elements with 

non-interior node F-27 

2. Plate problems with 

a) N £2 varying F-‘28, F-36 

b) N-| 1 varying N 22I 

1) Plate element J is affected 

1st, plate element is affected F-ll 

2) Plate element J is not affected 
1st, plate element is affected F-12 

3) Plate element J is affected 

1st, plate element is not 
affected F-13 

4) Plate element J is not affected 

1st, plate element is not 
affected F-14 



9.40 


Routine 


Purpose 


Configuration 


Tests 


UPRBND 


5) Biaxial ratio ZO, plate 
element J is affected, first 
plate element is affected 

6) Biaxial ratio ZO, plate 
element J is not affected, 
first plate element is 
affected. 

7) Biaxial ratio ZO, plate 
element 0 is affected, first 
plate element is not affected 

8) Biaxial ratio ZO, plate 
element 0 is not affected, 
first plate element is not 
affected 

UPRBND calculates the panel upper Problem with 


bound for LOADOP = 5, 6, 7, 8. 

1 . 

Start substructure only 


2. 

Start and end substructure only 


3. 

Start, repeat and end substructures 


4. 

LOADOP = 5 or 7 where 


a. Subelement is already buckled 

b. reduced subelement is already 
buckled 

c. workable 
LOADOP = 6 or 8 


F-15 


F-16 


F-17 


F-18 


B-l 

B-22 

B-21 

B-6 

B-7 

B-5 


5 . 



9.41 


Routine 


Purpose 


USERID Converts user ID's for nodes, plate 

elements and beam elements to 
internal ID's and vice versa. 

VIPDR VIPDR produces the inner product of 

2 vectors in double precision. 


Configuration Tests 

a. subelement f 0 and ^ i 0 B-9 

b. subelement = 0 and ^ f 0 B-10 

c. subelement f 0 and ^ = 0 B- 1 2 

d. subelement N-j-j = 0 and ^ ^ 0 B-ll 

e. element =0 and N 22 t 0 B-16 

f. subelement buckled B-13 

g. reduced subelement buckled B-14 


B-15 
B-l 7 
B-l 8 
B-l 9 
B-20 


1 . 

Invalid type for ID search 

Inspection 

2. 

Invalid search option 


3. 

User ID not found 

F-64 

4. 

Internal ID not valid 



F-l 



9.42 


Routine 

YB 

YY 


7J\RK 


Purpose Configuration Tests 

YB computes the particular solution Thermal problem with plate. S-l 

of a plate elemental nodal displacement 
at Y = -B/2 or Y = b/2 

YY generates inter-element forces Thermal problem with plate. S-l 

due to temperature change from 
particular solution for plate 
element along X-axis at 
Y = b/2 or Y = -b/2. 

To find N zeros of an arbitrary F-l 

complex - valued function of a 
complex variable. 



10.0 


PROGRAM VALIDATION 


The program, BUCLASP3 was validated by running various data sets such 
that all major logic paths were tested. The validation is divided 
into two categories: 

a. Functional Tests 

b. Inspection 

Listed below is the functional tests and their characteristics. 
Section S.3 states for each program/subroutine the test(s) used. 


10il 



10.1 

Input Option 

Functional Test Correlation 


Card 

Field 

Option 

Test 

1 

1 

Data Set Title 

F-l 

2 

1 

I PC ( 1 ) = 0 

F-l 



IPC(l) = 1 

F-81 


2 

I PC (2) = 0 

F-l 



I PC (2) = 1 

F-82 


3 

I PC (3) = 0 

F-l 



I PC (3) = 1 

F-83 


4 

IPC (4) = 0 

F-l 



I PC (4) = 1 

F-83 


5 

I PC ( 5 ) = 0 

F-l 



I PC ( 5 ) = 1 

F-84 


6 

I PC (6) = 0 

F-l 



I PC (6 ) = 1 

F-85 


7 

IPC (7) = 0 

F-l 



I PC (7) = 1 

F-86 


8 

I PC (8) = 0 

F-l 



I PC (8) = 1 

F-87 


9 

I PC (9) = 0 

B-l 



I PC (9 ) = 1 

B-23 


10*2 



Card 


Field 

Option 

Test 

10 

IPC(IO) = 0 

B-l 


IPC(IO) = 1 

B-23 

11 

IPC(ll) = 0 

B-l 


I PC ( 11 ) = 1 

B-23 

12 

I PC (12) = 0 

B-l 


I PC ( 1 2 ) = 1 

B-23 

13 

I PC (13) = 0 

B-l 


I PC (13) = 1 

B-23 

14 

I PC (14) = 0 

B-l 


I PC (1 4) = 1 

B-23 

1 

JPC(l) = 0 

F-l 


JPC(l) = 1 

F-41 

2 

JPC(2) = 0 

F-36 


JPC(2) = 1 

F-l 

3 

JPC(3) = 0 

F-l 


JPC (3) = 1 

F-41 

4 

JPC (4) = n 

F-89 

5 

J PC ( 4 ) = k 

F-37 


10.3 



Card 

Field 

Option 

Test 

3 

6 

Default thermal input = 1 

S-7 



Default thermal input = 2 

S-l 


7 

Default value for number of 




subelements 

B-23 


8 

First Fourier Harmonic for 
printi ng 

S-l 


9 

Last Fourier Harmonic for 
printing 

S-l 


10 

Output option for stress 
NUMX > 0 

S-l 



Output option for stress 
NUMX < 0 

S-8 

4 

1 

Number of Nodes 

F-l 


2 

Number of Beams 

F-41 


3 

Number of Flat Plates 

F-41 


4 

Number of Curved Plates 

F-41 


5 

Section Length 

F-l 


6 

Load or Strain Value 

F— 19 


7 

Biaxial Ratio 

F-15 


10.4 



Card 

Field 

Option 

Test 

5 

1 

Load Option = 4 

S-l 



= 5 

B-l 



= 6 

B-2 



= 7 

B-3 



= 8 

B-4 


2 

Wave Number Option = 1 

F— 37 



= 2 

F-l 



= 3 

F-36 



= 4 

F-39 


3 

Loop Start Value 

F-l 


4 

Loop End Value 

F-l 


5 

Lower limit for root search 

F-35 


6 

Upper limit for root search 

F-35 


7 

Default number of subdivision 




= blank 

F-l 



= n 

F-37 


8 

Lower bound load 

F— 37 


9 

Upper bound load 

F-37 

6 

1-16 

M-list values 

F-36 

7 

1-14 

M- list values 

F-37 



Field 

Option 

Test 

variable 

Thickness Table 

F-37 

variable 

Material Table 

F-37 

1-9 


F-27 

1-9 


F-27 

1 

Letter N 

F-2 

2 

User Node Number 

F-2 

3 

Nodal Boundary Condition 



= blank 

F-2 


= 1 

F-2 


= 2 

F-l 


= 3 

F-l 


= 4 

F-l 1 

4 

Y-coordinate 

F-2 

5 

Z-coordinate 

F-2 

1 

Letter "D" 

F-l 

2 

Number of substructures 

F-l 

3-8 

Substructure Node Pairs 



start only 

F-l 


start and end only 

F-38 


start, repeat, and end 

F-22 



Card 

S 


F 


Field 

Option 

Test 

1 

Letter "S" 

F- 11 

2 

User Node ID 

F-ll 

3 

w component 

F-ll 

4 

Q component 

F-ll 

5 

v component 

F-ll 

6 

u component 

F-ll 

1 

Letter "F" 

F-22 

2 

User Plate ID 

F-22 

3 

0FF1 

F-22 

4 

0FF2 

F-22 

5 

0FF3 

F-22 

6 

0FF4 

F-22 

1 

Letter "P" 

F-l 

2 

User Plate ID 

F-l 

3 

Node I 

F-l 


10.7 



Card 

P 


Field 

Option 

Test 

4 

Node J 

F-l 

5 

Plate Type = 1 Flat 

F-l 


= 2 Curved 

F-30 

6 

Element offset = 0 No 

F-l 


= 1 Yes 

F=22 

7 

Number of laminas 

F-l 

8 

Input mode = 0 P-1, or P-2 card 

F-l 


= 1 Table 

F-37 

9 

N 99 load switch = 0 affected 

F-38 


= 1 not affected 

F-36 

10 

Number of subdivisions 

F-37 

11 

Type of material properties 



= 0 Engineering Constants 

F-37 


= 1 Q matrix 

F-37 

12 

Local plate thermal input = 0 

S-l 


= 1 

S-7 


= 2 

B-23 


= 3 

B— 23 

13 

Load input option = 0 

B-3 


= 1 

B-24 

14 

Substructure assignment 

S-l 


10*8 



Card 

Field 

Option 

Test 

P 

15 

Number of subelements 

B-23 


16 

Curved Plate Radius 

F-21 

P-1 

1 

T 

F-l 


2 

E 11 

F-l 


3 

E 22 

F-l 


4 

RNUA 

F-l 


5 

G 1 2 

F-l 

P-2 

1 

T 

F-37 


2 

Qn 

F-37 


3 

Ql2 

F-37 


4 

Q 22 

F-37 


5 

Q 66 

F-37 

P-3 

1 

TX 

S-7 


2 

TYA 

S-7 


3 

TYB 

S-7 


4 

TYC 

S-7 


10.9 



Card 

Field 

Option 

Test 

P-3 

5 

TZA 

S-l 


6 

TZB 

S-l 


7 

TZC 

S-l 

P-4 

1 

ALPHX 

S-7 


2 

ALPHY 

S-7 


3 

TZ 

S-7 

P-5 

1-8 

i input 

B-l 

P-6 

1-8 

N^2 input 

B-l 

B 

1 

Letter "B" 

F-23 


2 

User Beam ID 

F- 23 


3 

Beam Node 

F-23 


4 

Type of Beam = 1 

F-25 



= 2 

F-24 



= 3 

F-26 


5 

Number of layers 

F-23 


10.10 



Field 

Option 

Test 

6 

Input Mode = 0 B-l , B-2 or 



B-3 card 

F-23 


= 1 By tables 

F-39 

7 

Thermal Input Option = 0 

S-5 


= 1 

S-9 

8 

Load Input Options = 0 

B-l 


= 1 

B-24 

9 

Beam Angle 

F-23 

10 

Beam Area 

F-23 

1 

E-Modulus 

F-25 

2 

G-Modulus 

F— 25 

3 

Moment about local yy-axis 

F-25 

4 

Moment about local zz-axis 

F-25 

5 

Warping constant 

F-25 

6 

Torsion constant 

F-25 

7 

Y offset 

F-25 

8 

Z offset 

F-25 


10.11 



Card 

Field 

Option 

Test 

B-2 

1 

Thickness 

F-23 


2 

E-Modulus 

F-23 


3 

G-Modulus 

F-23 


4 

Beam width 

F-23 

B-3 

1 

Radi us 

F-26 


2 

E-Modulus 

F-26 


3 

G-Modulus 

F-26 

B-4 

1 

MTY 

S-10 


2 

MTZ 

S-10 


3 

PIT 

S-10 

B-5 

1 

TXO 

S-ll 


2 

TAZ 

S-ll 


3 

TBZ 

S-ll 


4 

TCZ 

S-ll 


10.12 



Card 

Field 

Option 

Test 

B-6 

1 

ALPHX 

S-ll 


2 

TY 

S-ll 

B-7 

1 

ALPHX 

S-ll 


2 

TR 

S-ll 

B-8 

Axial load 


B-21 

L 

1 

Letter "L M 

F-l 


10.13 



10.2 


Functional Test Descriptions 


Test 


F-l 


F-2 


F-3 


F-4 


Characteristics 

1. Flat Panel at 90° SS - Free 

2. Buckling and eigenvector solution 

3. M0PT * 2 

4. Flat Plates with 

a. left node reduced out 

b. both nodes interior 

c. right node reduced out 

5. Input by P-1 cards 

6. Simply-supported node 

1. Flat Panel at 90° C - Free 

2. Clamped node 

3. Free node 

1. Flat Panel at 0° angle SS - Free 
1. Flat Panel at 0° angle C - Free 


F-5 

F-6 

F-7 

F-8 

F-9 


1. Panel at angle (non-orthogonal) SS - Free B.C. 
1, Panel at angle (non-orthogonal) C - Free B.C. 
1. Flat plate with non-zero B matrix 
1. Curved panel at angle (1) 

1. Curved panel at angle (2) 


F-ll 


1. varying 

2 . H22X l n P u * 


10.14 



10.2 F-ll Continued ... 


F-12 


F-13 


F-14 


F-15 


F-16 


F-17 


3. First Plate Affected 

4. J th Plate Affected 

5. Sprung B.C. 

(J is the critical element for panel upper bound.) 

1. N^ Varying 

2 . ^221 

3. First Plate Affected 

4. 0 th Plate Not Affected 

1. Varying 

2. Nggi Input 

3. First Plate Not AFfected 

4 . J th Plate Affected 

1. Varying 

2. N^2 j Input 

3. First Plate Not Affected 

4. J th Plate Not Affected 

1. Biaxial Ratio 

2. First Plate AFfected 

3. 0 th Plate Affected 

1. Biaxial Ratio 

2. First Plate Affected 

3. J th Plate Not Affected 

1. Biaxial Ratio 

2. First Plate Not Affected 

3. J th Plate Affected 


10*15 



]0.2 Continued ... 
Test 
F-18 


F-19 


F-20 


F-21 

F-22 

F-23 

F-24 

F-25 

F-26 


Characteristics 

1. Biaxial Ratio 

2. First Plate Not Affected 

3. J th Plate Not Affected 

1. N 22 Varying 

2. EPS 11 Input 

3. First Plate Affected 

4 . J th Plate Affected 

4 » ' 

1. N 22 Varying 

2. EPS 11 Input 

3. First Plate Not Affected 

4 . J th Plate Affected 

1. Mix of Curved and Flat Plates 

1. Offsets on Plates 

2. Substructure Interrelationship Element Pairs 

3. Plates with Different A, B & D Matrices 

1. Rectangular Beam with 2 Layers 

2. Angle of Beam is Non-Trivial 

1. Rectangular Beam with 1 Lamina 
1. General Beam Element 
1. Circular Beam with 3 Laminas 
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10.2 


Continued ... 


Test 

F-27 

F-28 

F-29 

F-30 

F-31 

F-32 


Characteristics 

1. Plate with 3 Layers 

2. Rectangular Beam with 15 Layers 

3. Similar Plates with Different Boundary Conditions 

4. Similar Plates with Same Boundary Conditions 
with Different Widths 

5. Non-Neglected Plate Element with Non-Interior 
Node 

1. Flat Plate with Zero B-Matrix 

2. N 22 Varying 

3. Non-Zero N^ Side Load Input 

1. Flat Plate with Zero B-Matrix 

2. Nj-j Varying 

1. Curved Plate with Non-Zero B-Matrix 

2. Non-Zero Side Load Input 

3. Ng 2 Varying 

1. Curved Plate with Non-Zero B-Matrix 

2. Non-Zero N 22 Side Load Input 

3. Nj. Varying 

1. Determinant Always Positive 

2. 0/2 Buckling Condition 


F-33 1.' Curved Elements with Zero B-Matrix and Different 

Radii 

2. Similar Plates with Same Boundary Conditions 
with Same Widths 

3. Similar Plates with Same Boundary Conditions with 
Same Widths and Same Spring Stiffness 
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10.2 


Continued 




Test 

F-34 

F-35 

F-36 


F-37 


F-38 


F-39 


F-40 


F-41 


Characteristics 

1. Curved Plate with Zero B-Matrlx, N 22 
Varying . 

1. Multiple Root Searching 

1. N 22 Varying, "11 Input Case with Some 

Unaffected Elements. M0PT * 3 Used and 
Buckling Only Option Used 

: - i, *• 

1, Flat Plate with Blj f 0 with Input Upper Bound 
and M0PT * 1. Qlj of One Plate Input with P-2 
Card, Qlj of Another Plate Input with Table and 
Card. E, v Input by Table, JPC(5) ■ 6; 

Default Subdivisions • 10; One Element Set * 5 

1. Two Identical Plates with Different Spring 
Stiffnesses at the Same Sprung Sides Going 
Through the Upper Bound Calculations. Only 
Start and End Substructures are Used for This 
Problem. 

1. One Beam Is Circular with One Layer, the Other 
Is a General Beam with Offsets. The Circular 
Beam Is Input with Table. M0PT ■ 4 Is Used. 

1. Two Laminated Plates with Same Loads and Radii, 
But with Different Widths. These Plates are 
Not Considered Identical From the Stiffness 
Standpoint. 

1, Too Many Nodes 
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10.2 


Continued 

Test 


i » i 


Characteristics 


F-42 

F-43 

F-44 

F-45 

F-46 

F-47 

F-48 

F-49 

F-50 

F-51 

F-52 

F-53 

F-54 

F-55 

F-56 

F-57 

F-58 

F-59 


1. Too Many Elements 
1. Too Many Beams 
1. Too Many Plates 

1. Bad Equivalent Node Data 

2. Equivalent Nodes with Different Boundary Conditions 
1. Bad Interrelationship Data 

1. Too Many Values in M-list 
1. No L Card 

1. Non-Element Data Out of Order 

1. Too Many Nodal Spring Stiffnesses 

1. Too Many Plate Offsets 

1. Too Many Interrelationship Element Pairs 

1. Too Many Equivalent Node Pairs 

1. Bad Table Input 

1. Incorrect Number of Nodes Input 

1. Bad Thickness Table 

1. Too Many P Cards 

1. Too Many B Cards Input 

1. Invalid Load Option 
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Continued 


• • • 


10.2 

Test 

F-60 

F-61 

F-62 

F-63 

F-64 

F-65 

F-66 

F-67 

F-68 

F-69 
F-70 
F-71 
F- 72 
F-73 
F-74 
F-75 


Characteristics 

1. Invalid Wave Number Search Option 

1 . Too Many N Cards 

1. Illegal Boundary Condition 

1. Wrong Number of Spring Stiffnesses Input 

1. Bad User ID on Nodal Stiffness Set 

' t '* ■ . 

t - t. •• ■ 

1. Bad User ID in Interrelationship Data 

1. Bad User ID in Equivalent Node Data 

V. Bad User ID In Substructure Definition 

1. Incorrect Number of Substructure Definition 
Node Pairs Input 

1. Substructure Definition Element Pair Out of Order 

1. Interrelationship Between Substructures Invalid 

1. Invalid Node Used on Plate 

1. Invalid Node on Beam Definition 

1. Incorrect Number of Plates Input 

1. Incorrect Number of Beams Input 

1. Node In Start Substructure with Attachment to 
Element to Repeat, But Not Part of a Connection 
Node Pair. 
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10.2 


Continued 


* * • 


Test 

F-76 

F-77 

F-78 

F-79 

F-80 

F-81 

F-82 

F-83 

F-84 

F-85 

F-86 

F-87 

F-88 

F-89 

F-90 


Characterist ics 

1. Substructures with No Interior Node 

1. Panel with 2 Substructure, But No End Sub- 
structure 

1. Boundary Condition on Beam Node 

1. Beam Node Which Is Not Equivalenced 

1. Element with Nodes In Start and End Sub- 
structure 

1. IPC(l) - 1 

!• I PC ( 2 ) - 1 

1. IPC(3) -1 

1. IPC(4) • 1 

1. IPC(5) - 1 r 

1. IPC(6) - 1 
1. IPC(7) - 1 
1, IPC(8) • 1 

1. JPC(4) - n < 100 and n f 0 

i 

1. Duplicate User 10 on Nodes 

2. ' Duplicate User ID on Beams 

3. Duplicate User ID on Plates 


10.21 



Thermal Stress Problems 


Characteristics 

1. Four flat plates 

2. Two general beams 

3. Simply supported node 

4. Plate with left node constrained 

5. Plate with right node constrained 

6. Plate with no node constrained 

7. Repeated element types 

8. Start substructure only 

1. Two curved plate plate with non-zero B 
matrices 

2. Repeated element types 

3. Sprung node 

See Section S-6 

1. Beam in repeat substructure 

2. Plate in repeat substructure 

1. Plate offsets 

2. Beam at angle 

3. Start, Repeat and End substructures 

1. Start and End substructures 

2. Curved plates with zero B-matrix 

3. Repeated element types 

4. Sprung node 

1. T z input 

2. Flat plate with non-zero B-matrix 



Test 


Characteristics 


S-8 NUMX < 0 

S-9 No thermal input on beam 

S-10 Thermal Data on General Beam 

S- 11 1. Thermal Data on Rectangular Beam 

2 . Thermal Data on Circular Beam 
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Test 


Characters istics 


B-l 


B-2 


B-3 


B-4 


B-5 


B-6 


B-7 


B-8 


B-9 


B-l 0 


B-l 1 


1. LOADOP = 5 

2. Two subelements per elements 
LOADOP = 6 

1. LOADOP = 7 

2. Start Substructure only 

LOADOP = 8 

Upper bound for LOADOP = 7 
N-ji f 0, N 22 t 0 and 

Subelement > Galerkin upper bound 

Subelement > Reduced subelement upper bound 

= 1000 constant 

Upper bound for LOADOP = 8 

All N-j ^ ^ 0 and N 22 f 0 

N-ji = 0 on one subelement, rest f 0 
N 22 = 0 on same subelement, rest + 0 

N-ji = 0 on one subelement, rest f 0 
"N 22 = 0 on same subelement, rest f 0 
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Test 


Characteristics 


B-12 

B-l 3 

B— 1 4 

B-l 5 

B-16 
B-l 7 
B-18 

B-19 

B-20 

B— 21 


i 


N 

N 


11 

22 


t 0 on one subelement, rest f 0 
= 0 on same subelement, rest f 0 


Subelement > Galerkin upper bound 

N ^2 = 0 for same subelement 

Subelement > Galerkin upper bound same 

Subelement > Galerkin upper bound] subelement 

Subelement > Galerkin upper bound 

N-j ^ = 0 for same subelement 


N.|.| = 0 an< * ^22 = ® f° r one e ^ ement 

Nii = 0 anc * ^22 = ® ^ or elements 

Subelement Nn > reduced subelement upper bound 
N ^2 “ 0 for same subelement 


Subelement Nn > reduced subelement upper bound 
Subelement N ^2 > reduced subelement upper bound 

same Subelement N^ > reduced subelement upper bound 

subelement Nn = 0 for same subelement 


1 . Beam Load 

2. Start, Repeat and End substructures 

3. - Externally loaded 
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Test Characteristic 

B-22 1. Externally loaded 

2. Start and End substructure 

B-23 1. I PC ( 9 ) = 1 

2. IPC(IO) =1 

3. IPC(ll) - 1 

4. I PC (12) = 1 

5. I PC ( 13) = 1 

6. I PC (14) = 1 

7. Input default number of subelement 

8. No thermal input on one plate 

B-24 No load on one plate element 
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APPENDIX 


During the running of the thermal buckling problem shown in Section 8.3, 
singularities were encountered in the program in the subelement processing 
stage. By re-modelling the cylinder, a successful solution was achieved. 

The original idealization consisted of three equal wide elements between 
nodes 8 and 21 of Figure 8.3. The re-modelled cylinder shown in Figure 8.3, 
has thirteen elements ( 8 to 20 ) instead of the original three elements. 

This change in idealization which reduces the element width, 'bypassed' the 
singularity problems. 



